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ABSTRACT 


\ 

Volume  Two  comprises  three  appendices  giving  detailed  analyses  of 
the  problem  areas  in  the  investigation  of  the  use  of  FDMA  with  a mix  of 
user  terminals  of  differing  characteristics.  Appendix  A deals  with  the 
limiter  transfer  characteristic  and  provides  a comprehensive  theoretical 
investigation  of  the  problem  of  transmitting  a number  of  constant-envelope, 
phase-modulated  FDMA  carriers  through  a limiting  satellite  repeater. 
Appendix  B analyzes  the  problem  of  detection  of  quadriphase-modulated  FDMA 
carriers  at  the  receiver  after  transmission  through  a limiting  transponder. 
An  expression  for  the  bit-error  probability  at  the  output  of  an  FDMA  chan- 
nel is  obtained  considering  the  presence  of  other  FDMA  carriers,  cross 
products,  and  retransmitted  up-link  as  well  as  down-link  noise  at  the  re- 
ceiver input.  Appendix  C describes  the  computer  program  SYSCON.  It  also 
discusses  the  equations  and  subroutines  used  in  system  optimization. 
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Appendix  A 


INVESTIGATION  OE  THE  EFFECTS  OF  BANDPASS  LIMITING 
ON  TRANSMISSION  OF  FDMA  SIGNALS 


1 . Introduction 

This  appendix  documents  the  results  of  the  investigation  performed 
to  provide  a comprehensive,  theoretical  understanding  of  the  transmission 
of  a number  of  constant-envelope,  phase-modulated  FDMA  carriers  through 
a limiting  satellite  repeater,  and  to  obtain  expressions  for  the  signals, 
intermodulation  products,  and  up-link  noise  at  the  repeater  output. 

The  model  generally  used  in  theoretical  investigations  to  represent 
a limiting  satellite  repeater  is  shown  in  Figure  A-l.  The  limiting  can 
be  either  hard  or  soft.  The  bandpass  filter  preceding  the  limiter  is 
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FIGURE  A-1  SATELLITE  REPEATER  MODEL  FOR  INVESTIGATION  OF 
IDEAL  SYMMETRIC  LIMITING 


wide  enough  to  pass  the  signals  with  negligible  distortion;  it  limits  the 
repeater  input  noise  to  a bandwidth  that  is  small  compared  to  the  center 
frequency  of  the  filter.  The  bandpass  filter  following  the  limiter  con- 
fines the  output  spectrum  to  essentially  only  the  fundamental  band  of  the 
signals. 
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The  model  of  1- inure  A- 1 requires  an  "ideal"  zonal  filter  at  the 
limiter  output,  which  passes  the  fundamental  band  of  signals  unattenuated 
but  completely  suppresses  the  higher  harmonics  oi  the  signal  frequencies. 
Heal  filters  can  only  approximate  this  desired  behavior.  The  assumption 
(in  Figure  A-l)  of  an  ideal  bandpass  filter  following  the  limiter  neg- 
lects principally  the  in-band  effects  of  the  filter  on  the  signals,  the 
intermodulation  products,  and  the  noise  present  at  the  limiter  output. 
However,  complete  analysis  of  a more  realistic  bandpass  filter  would  be 
extremely  difficult,  since  the  system  is  nonlinear  with  memory.  A 
rigorous  analysis  of  the  problem  would  involve  solving  lor  the  filter 
output,  which  is  related  to  the  limiter  input  process  by  a nonlinear  dif- 
ferential equation  determined  by  the  limiter  and  filter  characteristics. 

Section  2 of  this  appendix  develops  analytical  expressions  to  repre- 
sent the  nonlinear  transfer  characteristics  for  the  various  types  of 
limiters  (e.g. , quasi-1  inear , soft,  and  hard)  that  would  be  of  interest 
in  an  KDMA  system.  Section  3 provides  a comprehensive,  theoretical 
analysis  of  the  bandpass  limiter  model  of  figure  A-l.  Expressions  for 
determining  the  signal,  intermodulation,  and  noise  components  at  the  out- 
put of  the  satellite  repeater  are  developed  by  assuming  the  input  to  con- 
sist of  a sum  of  n signals  of  arbitrary  amplitude  and  phase  modulation. 
The  analysis  is  conducted  completely  in  the  time  domain,  to  retain  the 
phases  of  the  signals  and  intermodulation  products  at  the  limiter  output. 
This  is  not  possible  with  the  autocorrelation  function  approach  normally 
used  in  limiter  analysis.  Section  4 summarizes  the  results  of  Sections  2 
and  3 and  offers  conclusions. 
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Limiter  Character is  tics 


a.  General 

To  analyze  the  problem  of  transmitting  a number  of  constant- 
envelope,  phase-modulated  sinusoidal  carriers  through  a limiting  satellite 
repeater,  we  need  a mathematical  representation  lor  the  nonlinear  charac- 
teristic of  tne  limiter.  Generally,  three  types  of  limiter  characteris- 
tics are  of  interest  in  satellite  communications: 


• Hard- 1 imi t ing 

• Quasi-linear-limiting 

• Soft-limiting. 

These  characteristics  are  shown  in  Figure  A- 2.  In  the  hard-limiting  case, 

the  output  is  either  fry.  The  quasi-linear  limiter  has  a linear  region 
with  a sharp  break  at  the  saturation  point,  while  a soft  limiter  has  a 
gradual  transition  from  the  linear  region  to  the  saturation  region. 


b.  Mathematical  Representation 

The  mathematical  representation  for  each  of  the  three  limiter 
characteristics  in  Figure  A-2  may  be  obtained  from  an  infinite  integral 
that  relates  the  limiter  output  y to  its  input  x as 


0 

H 


sin(vx)  • sin(vtv/n)  e 


2 2 

-v  v dv 


( A- 1 ) 


where  ry  is  the  constant  output  limiting  level,  and  the  parameters  n and 
Y,  respectively,  determine  how  fast  saturation  is  approached  with  a 
quasi-linear-  and  a soft-limiting  characteristic. 

The  characteristic  of  a quasi-linear  limiter  is  expressed  in 
analytical  form  as 
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OUTPUT  V 


HARD-LIMITER 
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FIGURE  A-2  LIMITER  TRANSFER  CHARACTERISTICS 
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i-1 


s i n ( v x ) 


(IV 

s i n ( v n n ) — 


(A-2) 


which  is  obtained  from  Kq.  ( A— 1 ) by  sotting  v = l*.  Equation  (A-2)  is 
plotted  in  figure  A-2  and  shows  how  last  saturation  is  approached  as  n 
is  increased.  for  large  n,  the  quasi-linear  characteristic  becomes 
hard- 1 imi t ing. 


A soft  limiter  is  usually  described  by  an  error-function  charac 
teristic.  This  may  be  obtained  from  Eq.  (A-l),  if  n approaches  infinity 
so  that  sin(v<y  n)  vry  n,  and 
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OUTPUT  Y 


FIGURE  A-3  QUASI-LINEAR  LIMITER  VOLTAGE  TRANSFER  CHARACTERISTIC 


a 
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2 2 

2rv  | . , x -v  y /2  dv 

y = — I sin(vx)  e — 

n I v 


(A-3) 


Integration  of  Eq.  (A-3)  yields  the  following  more  conventional 
. * 

form  for  soft  limiting:1 


y = n erf ( \//2y) 


(A-4) 


References  are  listed  at  the  end  of  Volume  Two. 
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where  the  error  function  is  defined  as: 


( A-5 ) 


Equation  (A-4)  is  plotted  in  Figure  A-4  and  shows  how  saturation  is 
approached  as  the  parameter  y is  decreased.  In  the  limit,  when  y 
approaches  zero,  we  obtain  the  well-known  expression  for  the  character- 
istic of  a hard  limiter  [y  = 0,  in  Eq.  (A-3)]: 


c . 


o 


dv 

v 


Power  Output 


( A-6 ) 


In  practice,  the  different  limiter  characteristics  are  often 
described  in  terms  of  power  output  versus  power  input,  assuming  the  input 
to  be  a single  sinusoidal  carrier.  I t is  therefore  important  to  know  how 
the  power  output-input  characteristic  can  be  obtained  analytically  from 
the  voltage  output-input  characteristic  described  above  for  the  different 
limiter  types.  Below,  we  determine  the  output  power  for  the  three  types 
of  limiter  by  assuming  a sinusoidal  signal 


x = a cos  out 


(A- 7) 


at  the  limiter  input. 

By  substituting  the  above  expression  for  the  signal  in  Eq.  (A-l) 
and  using  the  Bessel  function  identity 


sin(z  cos  0)  = 2 'N  **  (-1)  • J (z)  cos(2m  + 1)  0 , (A-8) 

2m+l 


m=o 


the  following  expression  is  obtained  for  the  output: 
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y(t)  --n2J(-D 

m=o 


a 
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J ( va)  • sin(vry/n)  e 
2m+l 


2 2 

-v  V /2  dv 


cos(2m  + 1)  uut 


( A-9 ) 


The  output  at  the  fundamental  frequency  is  obtained  by  setting  m = o in 
the  above  expression: 


B(t) 


CO 

- n [ J (v 

" l 1 


2 2 

—V  v 2 dv 

a)  ■ sin(v»/n)  e — • cos  uat  . (A-10) 
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Expansion  of  the  sine  function  in  a power  series  yields: 


g(t) 


"E 


(-D 


p 


(2p  f 1).’ 


( va) 


2 2 

-v  y /2 


dv  cos  xt 


( A- 1 1 ) 


The  solution  of  the  above  integral,  which  is  attributed  to  Weber  and 
Sonine,2  is  given  in  terms  of  confluent  hypergeometr ic  series.  The  re- 
sult is: 


R<t> 


00  p 

2n  ^ (-1) 

Fry  0 / j ( 2p  + 1)  • pi 

p=o 


( A- 12 ) 


where  the  confluent  hypergeometric  series  is  defined  as 


a z a(  a + 1 ) 

iV3'  b,  z)  = 1 + - . - + — — 


2 

z 

2 : 


( A- 13 ) 


From  Eq.  (A-12)  we  can  obtain  the  output  of  each  of  the  three 
limiter  characteristics  and  can  also  calculate  the  power  output.  The 
output  of  a soft  limiter  is  obtained  by  letting  n approach  infinity. 
This  yields: 
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s(t)  = — a 
At 


COS  UQt 


(A- 14) 


where 


2 

P 


(A-15) 


The  power  of  the  sinusoidal  signal  at  the  input  is  given  by  P . The 

i 

power  output  of  a soft  limiter  is,  thus, 


P 

o 


2 

P 


,2/1 

1M2’ 


( A- 16 ) 


or 


2 - 

— 

2 2 2 -p 

P = “ a p e I 

/ 2 / \ 

1 2 / \ 

(p  /2 ) 

I,  (p  /2) 

O TT 

0 \ / 

1\  /J 

since 


/I  \ -z/2 

“ 

1F] 

lIS'  2’  -z)  = 6 

I (z/2)  f 
0 

I z/2) 

where  I ( ) and  I ( ) are  modified  Bessel  functions, 
o 1 

2 

For  large  values  of  p , the  confluent  hypergeometric  function 
simplifies  to3 


/I  2\  2 1 

lFl\2>  2>  ~P  ) ’ e**00  ’ (A_19) 

and  we  obtain  the  maximum  power  output  in  saturation,  which  is  the  same 
as  the  power  output  for  hard  limiting  (y  = 0): 


P 

max 


( A-20) 
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Equation  (A-20)  also  shows  the  relationship  between  the  hard- limiting 
power  output  and  the  limiting  level. 

Combining  Eqs.  (A- 16)  and  (A-20),  we  obtain  the  normalized 
power  output 


P 

o n 2 


max 


2/1 
1* i\2’ 


( A- 21) 


for  a soft  limiter  relative  to  that  for  a hard  limiter.  Equation  (A-21) 

can  also  be  viewed  as  representing  the  output  power  loss  resulting  from 

soft  limiting  compared  to  hard  limiting.  Equation  (A-21)  is  plotted  in 

2 

Figure  A-5  as  a function  of  p . 


The  output  of  a quasi-linear  limiter  is  obtained  by  letting  y 
approach  zero  in  Eq.  (A-12)  and  using  the  asymptotic  expansion  for  the 
continent  hypergeome tr ic  function:3 


g(  t)  = — a 


(-1) 


p=o  (2p  + 1) 


il  - p) 


,2p 

\ na  / 


cos  cut 


~S1  ( A-22) 

na 


where  F(  ) is  the  gamma  function.  The  series  in  Eq.  (A-22)  is 
as  can  be  easily  shown  by  d'Alembert's  ratio  test  for  absolute 
gence,  if  a/na  is  less  than  or  equal  to  one.  For  o/na  greater 
the  limiter  characteristic,  of  course,  is  linear  (Figure  A-3), 
have 


convergent , 
conver- 
than  one, 
and  we 


g(t)  = na  cos  uut  , — 5 1 

na 


The  output  power  of  a quasi-linear  limiter  is,  therefore, 


12 


2 2 
n a 

Po  = — 


for  — >.  1 
n a 


( A-23 ) 
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for  — £ 1 

na 


(A- 24) 


An  alternative  expression  for  the  quasi-linear  limiter  output, 
which  is  simpler  to  compute,  is  obtained  by  expanding  the  clipped  input 
waveform  in  a Fourier  series.4  The  output  signal  component  is  given  by: 


g(t)  = 


— — \ /l  - 4 p arc  sin^— j cos  tut  , for  p S 1 , (A-25) 

L * p p _ 


22  o p 

2 n a 2 i 

p = = 2n  — 

Ho  O 


Combining  Eqs.  (A-25)  and  (A-20),  we  obtain  the  normalized  power  output: 


P 4 

max 


1 1 2 + P UrC  Si 


ln(;) 


p ^ 1 (A- 26) 


for  a quasi-linear  limiter  relative  to  that  of  a hard  limiter.  Equation 

2 

(A- 26)  is  plotted  in  Figure  A-6  as  a function  of  p . It  is  interesting 

to  note  that  when  the  amplitude  of  the  input  signal  a is  equal  to  <y  n, 

2 

i.e.,  p = 1,  the  output  power  is  smaller  by  a factor  of  tt  16  (-2.1  dH) 
than  the  maximum  power  obtained  at  saturation. 
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Practical  Limiter  Characteristics 


Hard  or  quasi-linear  limiting  is  rarely  encountered  in  practice. 
A soft  limiter  with  gradual  saturation  intuitively  appears  to  most  closely 
characterize  a practical  limiter.  Generally,  the  power  output  of  prac- 
tical limiters  increases  linearly  (unity  slope)  at  low  input  power  levels 
and  reaches  saturation  at  large  power  levels.  figure  A-7  shows  the 
measured  characteristic  of  a tunnel  diode  limiter.  It  is  apparent  that, 
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FIGURE  A-7  MEASURED  POWER  TRANSFER  CHARACTER- 
ISTIC OF  A TUNNEL-DIODE  LIMITER 

to  approximate  closely  the  transfer  characteristics  of  practical  limiters, 
the  analytical  representation  should  provide  unity  slope  at  low  signal 
levels,  with  an  additional  parameter  available  to  characterize  the  soft- 
ness with  which  the  limiter  reaches  saturation. 
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The  slope  of  the  generalized  limiter  characteristic  of  Eq. 


( A- 1 ) is: 


a 

— = — n f 

dX  TT  J 


cos(vx)  • sin(vcy/n)  e 


2 2 

-v  v /2  dv 


Integration  of  the  expression  yields: 


dy  n I n 
— = — erf  I — z=— 
dx  2 I HZ\( 


at  low  signal  levels  (x  — o), 


- x 

n In 
— e r f I — — — 

2 1 /2y 


(A-27) 


<ly  a \ 

— = n erl(-rr — 1 
dx  \/ZnY/ 


( A-28 ) 


The  limiting  level  <y  is  a fixed  parameter  determined  by  the  specified 
power  output  when  saturation  is  reached.  From  Eq.  (A- 20),  it  is 


2 TT 
cy  = — P 


8 max 


(A- 29) 


Requiring  unity  slope  at  low  signals  yields 


" er,(/i^)  ■= 1 


( A-30) 


Since  the  error  function  can  never  exceed  one,  n must  be  one  or  greater. 
Consider  the  two  limiting  cases  for  n.  If  y = 0,  n must  be  equal  to  one, 
since  erf(oo)  = 1.  This  represents  the  transfer  characteristic  of  a 
quasi-linear  limiter  with  unity  slope.  Figure  A-8  shows  the  output/input 
power  characteristic  of  a quasi-linear  limiter  with  unity  slope  (n  = 1) 
at  low  power  levels.  The  limiting  level  a is  chosen  to  provide  0-dI3  out- 
put power  at  saturation.  If  n -*  oo,  the  argument  of  the  error  function 
will  tend  to  zero.  Expanding  the  error  function  in  a power  series  and 
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FIGURE  A-8  POWER  TRANSFER  CHARACTERISTICS  OF  QUASI-LINEAR,  SOFT.  AND  TUNNEL-DIODE  LIMITERS 


retaining  only  the  first  term  yields: 


Equation  (A-31)  is  the  relationship  for  determining  Y for  a soft  limiter, 
if  unity  slope  is  required  at  low  signal  levels.  figure  A- 8 also  shows 
the  output  input  power  characteristic  of  a soft  limiter  with  unity  slope 
at  low  input  power  levels. 

The  first  two  curves  in  Figure  A-8  correspond  to  the  limiting 
values  of  n (1  and  co).  lor  all  intermediate  values  of  n,  y can  be  de- 
termined from  Kq.  (A-30).  The  curves  representing  the  output/input  power 
characteristic  can  be  obtained  from  Eq.  (A-12);  they  will  all  lie  between 
these  two  curves. 

Comparison  of  the  soft-limiter  curve  in  Figure  A-8  with  the 
curve  in  Figure  A-7  shows  that  the  tunnel  diode  limiter  has  a very  soft 
saturation  characteristic  which  cannot  be  approximated  by  Eq.  (A-l). 

At  O-dB  input  power  level,  the  output  power  of  the  tunnel  diode  limiter 
is  approximately  -4.0  dB,  while  a soft  limiter  having,an  error-function 
characteristic  yields  -2.2  dB. 

Thus,  if  the  limiting  is  very  soft,  representation  of  the 
limiter  characteristic  in  terms  of  one  error  function  |Eq.  ( A—  4 ) | is  not 
sufficient.  Investigation  showed,  however,  that  the  use  of  two  additive 
error  functions  results  in  sufficient  flexibility  to  describe  even  very 
soft  limiter  characteristics.  The  limiter  characteristic  is  expressed 
as : 


(A- 3 2) 
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The  parameters  \ and  v v can  be  appropriately  chosen  to  approximate  the 
measured  characteristic  closely.  To  determine  y and  we  need  two 

independent  equations.  These  can  be  obtained  as  follows.  lor  small 
values  of  \,  differentiation  of  Eq.  (A-32)  yields: 


dy 

lx 


for  x — 0 


( A-33) 


Requiring  unity  slope  at  low  signal  levels  provides  the  first  condition: 


J 

Y 


1 


(A-34) 


A second  equation  is  obtained  by  first  specifying  an  output  power  for  a 
given  input  power.  lor  a sinusoidal  input  x = a cos  u)t,  the  output  is 
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( A- 3 5 ) 


where 


Pa 


2 P 

a i 

2y  2 = ~ 
2 v 
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( A- 36 ) 


Equation  (A-35)  is  obtained  by  summing  the  output  of  two  error  functions; 
the  output  of  one  error  function  is  given  by  Eq.  (A-14).  The  output  power 
is  t hen : 


P 

o 


( A- 3 7 ) 


Specification  of  an  output  power  P 

ol 

second  required  equation: 


at  an  input  power  P yields 
io 


the 
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( A- 38 ) 


Now,  and  Y,,  can  be  determined  by  solving  Eqs.  (A-34)  and  (A-38).  Ex- 
plicit expressions  for  Y^  and  Y,,  are  not  possible  because  of  the  confluent 
hypergeometric  functions  involved  in  Eq.  (A-38).  However,  they  can  be 
determined  numerically  by  solving  the  two  equations  on  a computer. 

The  curve  labeled  TDL  (tunnel  diode  limiter)  in  Figure  A-8  shows 
the  power  transfer  characteristic  based  on  Eq.  (A-37).  The  parameters  Y^ 
and  Yq  were  determined  by  specifying  an  output  power  level  of  -4.0  dB  at 
0-dB  input  power  level,  as  shown  by  the  measured  TDL  characteristic  in 
Figure  A-7.  Comparison  of  the  TDL  curves  in  Figures  A-7  and  A-8  shows 
excellent  agreement  between  the  measured  and  theoretical  results. 


e.  Summary 

The  voltage  transfer  characteristics  of  a quasi-linear,  a soft, 
and  a hard  limiter  can  be  expressed  analytically  in  terms  of  infinite 
integrals  by  relating  the  limiter  output  to  its  input  as  summarized  in 
Table  A-l.  In  the  case  of  soft  limiting,  the  integral  representation  can 
also  be  expressed  in  terms  of  an  error  function.  The  power  transfer 
characteristics  of  these  limiter  types  are  also  shown  in  Table  A-l. 

In  practice,  the  power  output  of  limiters  generally  increases 
linearly  with  unity  slope  at  low  input  power  levels  and  reaches  satura- 
tion at  high  power  levels.  The  use  of  a quasi-linear  or  a soft  limiter 
is  restricted  to  producing  a maximum  of  -2.2-dB  power  back-off  at  0-dB 
input  power  consistent  with  unity  slope  at  small  input  power  levels  and 
0-dB  output  power  at  saturation.  To  achieve  sufficient  flexibility  to 
describe  limiters  with  substantial  back-of f — e. g. , more  than  -2.2  dB  at 
0-dB  input  power — a slightly  different  representation  is  used  for  the 
limiter  characteristic.  The  characteristic  of  such  limiters  can  be 
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ANALYTICAL  REPRESENTATION 


Limiter  output  power  level  in  saturation 


represented  by  the  sum  of  two  soft  limiters  each  having  an  error-function 
characteristic  as  shown  in  Table  A-l. 

3.  Bandpass  Limiter  Analysis 
a.  General 

The  most  commonly  used  approach  for  determining  the  limiter 
output  components,  generally  known  as  the  characteristic  function  method, 
involves  computing  the  autocorrelation  function  of  the  limiter  output 
and  then  taking  the  Fourier  transform  to  obtain  the  power  density  spec- 
trum.6 Although  a general  expression  for  the  limiter  output  autocorrela- 
tion function  can  be  derived,  its  computation  becomes  extremely  involved 
when  modulation  of  the  input  signals  is  considered.  The  difficulty  lies 
in  determining  the  characteristic  function  of  the  signals  with  arbitrary 
modulation.  However,  if  the  signals  are  statistically  independent, 
angle-modulated  sinusoids,  and  only  the  average  power  or  the  magnitude 
of  the  signal  and  cross-product  terms  at  the  limiter  output  is  of  interest, 
the  modulation  of  the  input  signals  can  be  ignored,  which  considerably 
simplifies  the  analysis.  Davenport6  was  first  to  use  this  approach  to 
investigate  the  effect  of  hard  limiting  a single  unmodulated  sinusoidal 
signal  and  narrow-band  Gaussian  noise.  Jones  used  the  same  method  to 
analyze  the  case  of  two  sinusoids  plus  noise.  More  recently,  Shaft8  and 
Gyi3  independently  extended  the  analyses  to  include  n sinusoids. 

The  autocorrelation  function  method  has  the  inherent  drawback 
that  the  phase  information  of  the  signals  and  intermodulation  components 
is  lost  when  the  autocorrelation  function  of  the  limiter  output  is  de- 
termined. This  approach  is  therefore  not  convenient  for  FDMA  analysis, 
since  the  phases  of  the  limiter  output  components  are  also  of  interest. 

For  example,  in  the  computation  of  the  signal-to-intermodulntion-powor 
ratio  for  a given  channel,  the  phases  of  the  intermodulation  components 
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falling  in  that  channel  are  required,  since  the  spreading  of  the  power 
spectral  density  of  any  one  cross-product  term  is  a function  of  its  phase 
modulation,  even  though  the  average  power  in  that  output  component  is 
independent  of  its  phase.  Another  serious  disadvantage  of  using  the 
au tocorrela t iott  approach  is  that,  since  it  yields  only  the  autocorrela- 
tion function  of  the  limiter  output,  the  amplitude  distribution  of  the 
interference  (up-link  noise  and  cross  products)  cannot  be  determined. 

To  determine  the  bit-error  probability  in  the  detection  of  digitally 
phase-modulated  signals,  it  is  important  to  know  the  distribution  of  the 
noise  and  cross-product  terms  at  the  limiter  output. 

Stanford  Research  Institute  has  recently  developed  a time-domain 
method  of  analysis'  that  differs  markedly  from  the  older  mathematical 
techniques  widely  used  in  the  analysis  of  limiters.  The  main  feature  of 
this  analytical  approach  is  that  it  permits  the  limiter  output — consisting 
of  signals,  intermodulation  products,  and  noiso — to  be  expressed  com- 
pletely in  the  time  domain.  Representation  in  the  time  domain  permits 
retention  of  the  phase  of  the  FOMA  signals  and  intermodulation  products 
at  the  limiter  output.  The  time-domain  technique  also  permits  evaluation 
of  the  distribution  of  the  interfering  signals  (cross  products  and  noise) 
at  the  limiter  output.  We  use  this  analysis  approach  for  determining  the 
bandpass  limiter  output  components  and  their  distribution. 

b.  Calculation  of  the  Limiter  Output 

The  specific  model  for  the  bandpass  limiter  to  be  considered 
has  been  shown  in  Figure  A-l.  The  input  to  the  limiter, 

x( i)  = s ( t ) n(t)  , ( A-39 ) 

consists  of  the  signal  s(t)  and  a band  of  zero-mean  stationary  Gaussian 
noise  n(t).  It  is  assumed  that  the  bandpass  filter  preceding  the  limiter 
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is  wide  enough  to  pass  the  signal  with  negligible  distortion 
the  input  noise  to  a narrow  bandwidth  that  is  small  compared 
frequency  of  the  filter.  The  limiter  is  followed  by  another 
filter  that  confines  the  output  spectrum  essentially  only  to 
mental  band  of  the  signal. 


and  limits 
to  the  center 
bandpass 
the  funda- 


It  is  assumed  that  the  limiter  has  a soft-limiting  character- 
istic, described  by  an  error-function  relationship.  Thus,  if  the  limiter 
input  is  x(t),  the  output  y( t)  may  be  expressed  in  analytical  form  by 
Eq.  (A-3)  as 


a 

I 


y(t)  = — a I sin|vx(t)|  e 

Tt 


2 2 

■v  y 2 dv 


(A- 40) 


The  limiting  level  a is  related  to  the  maximum  power  output  of  the  re- 
peater by  Eq.  (A-20) 


a 


ts-^ 


( A-  '4 1 ) 


Here  we  analyze  the  problem  with  a soft  limiter;  the  same  approach, 
however,  is  applicable  to  the  other  limiter  characteristics  considered 
in  Section  2.  We  therefore  provide  the  required  expressions  for  the 
signal,  intermodulation,  and  noise  components  at  the  output  of  the  band- 
pass limiter  for  the  other  limiter  characteristics  without  derivation. 
When  Expression  (A-39)  is  inserted  into  Eq.  (A-40),  the  sine  of  a sum  may 
be  expanded  into  the  sum  of  two  products  of  sine  and  cosine;  thus, 
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The  narrow-band  Gaussian  noise  at  the  limiter  input  may  be  expressed  as 


All  the  terms  in  the  above  equation  are  random  functions  be- 
cause of  the  presence  of  the  noise  envelope  r and  the  phase  <p.  However, 


only  the  first  integral  will  yield  an  average  output,  while  the  two  other 
integrals  will  not  contribute,  since  all  the  terms  contain  the  random 
phase  of  the  noise. 

The  average  limiter  output  at  any  arbitrary  time  t is  obtained 
by  averaging  over  all  possible  noise  amplitudes: 


z(t)  = E [ y ( t ) 


CO 

~ n j s i n [ v s(t)]  • E^J^vr)] 


2 2 

-v  y /2  dv 


(A-46 ) 


Since  the  noise  envelope  has  a Rayleigh  distribution 


P(r>  = e*p(-  » 

a \ 2a  / 


( A-47 ) 


where  a is  the  total  noise  power  at  the  limiter  input,  the  expected 

value  of  J (vr)  is  given  by 
o 


r -i  f / vV\ 

:(JJo<vr>j  = I JQ(vr)p(r)  dr  = expl-  — - — I 


(A-48) 


Substitution  of  Eq.  (A-48)  in  Eq.  (A-46)  yields 


z(  t ) = - 
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( A-49 ) 


The  above  integral  may  also  be  written  in  terms  of  the  error  function1  as: 


z(  t)  = 


s(t) 


( A- 50 ) 


The  corresponding  expressions  for  the  case  of  a hard  limiter  are  obtained 
by  setting  y = 0 in  Eqs.  (A-49)  and  (A-50).  The  presence  of  noise  a1,  the 
input  of  a soft  limiter  thus  tends  to  make  the  limiting  even  softer, 
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because  of  the  gradual  saturation  characteristics  of  the  error  function. 
This  was  also  observed  by  Jones.'  The  effect  of  soft  limiting  on  the 
average  output  is  mathematically  identical  to  adding  noise  to  a hard 
limiter,  provided  that  an  error-function  representation  can  be  used  for 
the  soft- limiter  characteristic. 

The  components  that  collectively  make  up  the  limiter  output 
noise  spectrum  are  the  second  and  third  integrals  of  Eq.  (A-45): 
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( A-5 1 ) 


Expressions  UK-49)  and  (A-51)  are  the  two  fundamental  equations 
for  the  analysis  of  the  limiter  in  the  time  domain.  In  Davenport’s 
notation,6  Eq.  (A-49)  represents  the  output  signal  and  intermodulation 
components  resulting  from  the  interaction  of  the  signal  with  itself 
(s  x s terms).  Similarly,  Eq.  (A-51)  represents  output  noise  components 
that  are  due  to  the  interaction  of  the  input  noise  with  itself  (n  X n 
terms)  and  with  the  input  signal  (s  x n terms).  Separation  of  the  noise 
output  into  n x n terms  and  s y n terms  is  not  possible,  without  knowledge 
of  the  signal  waveform  s(t). 


If  the  characteristic  function  method  had  been  used  to  calcu- 
late the  limiter  output,  the  phases  of  all  the  components  at  the  limiter 
output  would  have  been  lost,  and  a separation  of  the  output  components 
as  in  Eqs.  (A-49)  and  (A-51)  would  have  been  possible  only  in  terms  of 
the  autocorrelation  function  of  the  individual  components. 
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c.  Limiting  of  n Signals  anti  Noise 
1)  General 

The  signal  and  intermodulation  products  at  the  limiter 
output  are  obtained  by  inserting  the  expression  lor  the  limiter  input 
signals 


s(t)  = a.  cos[a).t  4-  ® ( t )j  , 


( A-52 ) 


i =1 


where  a^  and  o^(t)  represent  the  amplitude  and  phase  modulation  of  the 
FDMA  carriers,  into  Eq.  (A-49),  and  using  the  general  relationship: 
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(A-53) 


where  the  symbol  77  indicates  that  all  n of  the  Jn  coefficients  are  mul- 

Hi 

tiplied  together.  Thus,  the  limiter  output,  prior  to  bandpass  filtering, 
is  given  by 


z(  t) 
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( A- 55 ) 


represents  the  amplitude  coefficient  of  the  output  signal  and  intermodu- 
lation components. 

Note  that  the  modulation  of  the  signals  affects  only  the 
phase  of  the  intermodulation  components,  as  can  be  seen  from  Eq.  (A-54), 
and  that  the  amplitudes  of  the  signals  and  intermodulation  components 
|Eq.  (A-55)J  are  independent  of  modulation.  Equation  (A-55)  is  dependent 
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only  on  the  characteristic  of  the  limiter  and  the  amplitudes  of  the  sig- 
nals at  the  limiter  input.  Such  a separation  between  the  amplitude  and 
phase  of  the  limiter  output  components  is  extremely  convenient,  particu- 
larly for  the  analysis  and  computation  of  the  intermodulation  components 
The  significance  of  the  parameters  in  Eq.  (A-55)  is  as  follows: 


• n is  the  total  number  ol  signals  (accesses). 

• The  value  of  p is  the  harmonic  ol'  the  i**'  signal 

i 

appearing  in  the  cross  product  under  considera- 
tion. Thus,  p =0  indicates  that  the  i**'  signal 
does  not  contribute  to  the  cross  product,  while 
p =1  indicates  that  the  fundamental  appears, 
and  so  forth.  The  composition  of  the  cross 

product  is  indicated  by  the  subscript  p p ...p  . 

1 2 n 

A subscript  100... 0 indicates  that  the  output 
term  under  consideration  is  Signal  1.  A sub- 
script 2100... 0 indicates  that  the  second  har- 
monic of  Signal  1 is  mixed  with  the  fundamental 
of  Signal  2.  and  so  forth. 


• a is  the  amplitude  (voltage)  of  the  i input 

1 

signal. 

2 

• a is  the  noise  power  at  the  input  to  the  limiter. 

2 2 th 

Thus,  a ' 2a  is  the  SNR  of  the  i signal  at  the 
i 

input  to  the  limiter. 


2)  Amplitude  of  I.imiter  Output  Components 

The  limiter  output,  given  by  Eq.  (A-54)  , was  derived  by 
assuming  soft  limiting.  In  the  case  of  hard- 1 imi t ing  and  quasi-linear- 
limiting  characteristics  (Figure  A-3),  only  Eq.  (A-55)  will  be  modified. 


I 


The  expression  for  the  case  of  a hard  limiter  is  obtained  by  setting 
V = 0 in  Eq.  (A-55).  Eor  the  case  of  a quasi-Xinear  limiter  with  unity 
slope  (n  = 1 and  y =0),  the  following  expression  for  the  amplitude  co- 
efficients is  obtained: 


. (A-56) 


When  the  limiting  is  very  soft,  the  limiter  characteristic 
can  be  described  by  a sum  of  two  soft  limiters,  represented  by  error- 
function  characteristics  with  different  slope  parameters.  Thus,  if  the 
limiter  characteristic  is  represented  by  Eq.  (A-32),  the  limiter  output, 
prior  to  bandpass  filtering,  will  be  given  by  Eq.  (A-54).  The  amplitude 
coefficient  in  this  case  is  equal  to  the  sum  of  the  amplitude  coefficients 
of  two  error  functions  of  slope  parameters  and  y ^ , respectively, 


( A- 57 ) 


A closed-form  solution  for  the  amplitudes  of  the  output 
components  is  extremely  difficult  for  n larger  than  three,  Eor  n equal 
to  three  or  less,  the  amplitudes  of  the  signals  and  intermodulation  com- 
ponents can  be  obtained  in  terms  of  hypergeometric  functions  by  using  the 
expressions  derived  in  Ref.  10.  A much  more  straightforward  approach 
that  is  valid  for  all  values  of  n would  be  to  determine  the  amplitudes  by 
numerically  evaluating  the  integral  of  Eqs.  (A-55)  and  (A-56).  A computer 
program  for  evaluating  Eqs.  (A-55)  and  (A-57)  is  discussed  in  Appendix  t'. 
This  program  can  easily  be  extended  to  compute  Eq.  (A-56)  by  including 
the  additional  sine  term. 
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3)  Frequency  and  Phase  of  Limiter  Output  Components 

Equation  (A-51)  shows  that  the  limiter  output  will  contain 
components  whose  frequency  and  phase  are  given  by: 


f = P f.  Pi,  ...  Pi' 
11  2 2 n n 


0(t)=p0(t)  p c ( t ) ...  I p © ( t ) i—  |P  p 

11  2 2 n n 2 \ 1 2 


P ) . (A-58) 

n / 


where 

P,  , P9,  • • • > P (0,1,  2,...). 

12  n 

Obviously,  these  relationships  will  remain  the  same  whether 
hard  or  soft  limiting  is  used;  in  fact,  they  hold  for  all  odd-power  de- 
vices. We  make  the  following  definitions: 

• The  zone  of  an  intermodulation  product  is  the  sum 
of  the  p coefficients: 


Zone  = p F p F ...  p 
12  n 


( A- 59 ) 


• The  order  of  an  intermodulation  product  is  the 

sum  of  the  magnitudes  of  the  p coefficients: 

i 

Order  = ! p I + I p I + . . . + I p I . (A-60) 

1 1 2 1 1 n 1 

Here,  we  are  interested  only  in  the  terms  that  fall  in  the 

first  zone, 


( A-61) 
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since  all  other  terms  will  be  filtered  out  by  the  bandpass  filter  fol- 
lowing the  limiter. 


It  is  apparent  that  the  determination  of  all  the  terms 

that  will  fall  in  the  pass  band  of  the  filter  is  quite  a formidable  task 

when  n is  large.  Fortunately,  Shaft"  has  shown  that  the  dominant  output 

components  are  the  signal  and  cross-product  terms  for  which  Ip. I is 

max 

small  (-2).  This  approximation  considerably  simplifies  Kq.  (A-54), 

since,  instead  of  considering  all  the  p from  -oo  to  x,  it  is  necessary 

i 

to  consider  them  only  in  the  range  from  -2  to  2.  The  significant  filter 
output  components  may  be  determined  by  first  considering  the  different 
frequency  combinations  in  Eq.  (A-54)  and  then  selecting  the  desired  ones, 
i.e.,  the  combinations  that  will  yield  frequencies  in  the  pass  band. 

Since  the  limiter  lias  an  odd-order  transfer  characteristic, 
only  the  third-  and  fifth-order  intermodulation  components  are  i lie  signifi- 
cant contributors  to  the  output  distortion.  As  can  be  seen  from  Eq.  (A-60), 
for  the  third-order  intermodulation  products,  the  sum  of  the  magnitudes 
of  the  p coefficients  is  three,  while  for  the  fifth-order  cross  products 
it  is  five.  The  frequencies  of  the  third-  and  fifth-order  cross  products 
are  shown  in  Table  A-2. 

4 ) Determination  of  Limiter  Output  Noise  Components 

The  components  that  collectively  make  up  the  soft-limiter 
output  noise  may  be  obtained  by  substituting  Eq.  (A-52)  for  s(t)  in  Eq. 
(A-51),  and  expanding  sin[v  s(t)|  and  cos[v  s(t)|  with  the  aid  of  Eq. 

(A-53).  The  expansion  of  cos  fv  s(t)l  will  contain  cosine  terms  instead 
of  the  sine  in  Eq.  (A-53).  The  resulting  expression  then  represents  the 
combined  output  noise  produced  because  of  the  interaction  of  the  input 
noise  with  itself  (n  X n terms)  and  that  produced  as  a result  of  the 
interaction  of  the  input  signal  with  noise  (s  X n terms).  In  this 
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Table  A- 2 


FREQUENCIES  OF  THIRD-  AND  FIFTH-ORDER  CROSS  PRODUCTS 


Order 

Frequency 

3 

f 

f - 

f = f 

a 

b 

c 111 

3 

2 f 

- f 

= f 

a 

b 

21 

5 

f 

f 

f - f - f = f 

a 

b 

c d e 1 1 1 1 1 

5 

2f 

f 

- f - f = f 

a 

b 

c d 2111 

5 

f 

f + 

1 

II 

a 

b 

c d 1112 

5 

2f 

f 

- 2f  = f 

a 

b 

c 212 

5 

3f 

- f 

- f = f 

a 

b 

c 311 

5 

3f 

- 2f 

= f 

a 

b 

32 

expression,  the  term  corresponding  to  m and  all  p.  equal  to  zero  repre- 
sents the  direct  feedthrough  noise  (n  X n terms),  and  all  the  remaining 
terms  constitute  the  noise  produced  by  the  mixing  of  the  signal  with  the 
limiter  input  noise  (s  x n terms).  The  expression  for  the  direct  feed- 
through noise  (n  x n terms)  is  given  by 


4 

CO 

r 

J (vr) 

TT j (va ) 

2 2 

-v  Y /2  dv  r 
e — cos  jj  t 

> cp(  t rj 

. ( A-62) 
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Davenport0  has  shown,  for  the  case  ol  one  signal  and  noise,  that  only 


the  n X n terms  contribute  significantly  to  the  output  noise  at  low 


limiter  input  SNRs.  When  the  number  of  signals  is  sufficiently  large 


that  the  amplitude  distribution  of  their  sum  is  approximately  Gaussian, 


the  SNR  for  any  one  signal  will  be  small  at  the  limiter  input.  Conse- 


quently, from  Davenport's  results,  it  is  reasonable  to  expect  that  only 


the  direct  feedthrough  noise  will  contribute  significantly  to  the  limiter 


output  noise. 


By  using  an  exponential  approximation  for  the  product  of 


the  Bessel  functions 


TKK)  --H”)  < 


(A- 63) 


( A-64 ) 


represents  the  total  power  of  all  the  input  signals,  Eq,  (A-62)  can  be 


readily  integrated  to  yield: 


n <t>  . ilR- 

nXn  \J  rr  a 


1 r ( t ) r~  ~i 

ii’i  -,  2,  - — cos^t  ■ cp(t)J  , ( A-65 ) 


2 2 2 


where  a = Y • ' . 

3 2 


In  FDMA,  the  total  power  in  all  the  input  signals  is  much 


larger  than  the  power  in  the  up-link  noise  at  the  limiter  input.  The 


confluent  hypergeome tr ic  function  in  Eq.  (A-65)  may  be  set  equal  to  unity 


to  a very  good  approximation.  This  can  be  shown  by  calculating  the  ex- 


pectation of  the  hypergeometric  function: 


] 


■ 


where  p(r)  is  given  by  a Rayleigh  distribution  |Eq.  (A-47)],  and  a is 
the  total  noise  power  at  the  limiter  input.  Because  in  1'DMA  ■ is  very 
much  less  than  one,  the  series  in  Eq.  (A-66)  can  be  approximated  by  the 
first  term,  yielding  a value  of  unity  for  the  expectation  of  the  confluent 
hypergeome tr ic  function.  With  this  approximation,  Eq.  (A-65)  simplifies 
to : 


(A-68) 


Equation  (A-68)  shows  that  the  bandpass  limiter  output  noise  will  tend  to 
have  a Gaussian  distribution  of  variance 


It  can  be  seen  that  Eq.  (A-68)  represents  the  dominant 
output  noise  in  the  time  domain.  This  representation  is  particularly 
useful  for  the  evaluation  of  the  system  following  the  bandpass  limiter. 
Equation  (A-68)  may  also  be  represented  in  the  form 


T1  (t)  = x (t)  cos  U)t-y(t)  sin  jj  t 
n*n  1 o 1 o 
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(A- 70) 


where  x^(t)  and  y^(t),  the  amplitudes  ol  the  "in-phase  and  "quadrature' 
components,  are  given  by: 


[2  I a \ 

x^t)  = ~ l-*  ) r(t^  cos 


y^t)  = Y n 1 r(t)  Sin  r-^l) 


( A- 7 1 ) 


(A-72) 


The  functions  x^(t)  and  y^(t)  are  Gaussian  random  variables  having  zero 
mean  and  equal  variance  given  by  Eq.  (A-69). 

Expressions  (A-68)  through  (A-72)  are  also  applicable  for 

2 

a hard  limiter  (y  =0).  The  only  change  will  be  to  replace  cs~  in  these 
2 

expressions  by  If  the  limiting  is  very  soft,  so  that  the  limiter 

characteristic  is  described  by  the  sum  ol  two  soft  limiters  represented 
by  error-function  characteristics  with  different  slope  parameters,  the 
bandpass  limiter  output  noise  can  be  obtained  by  summing  the  noise  con- 
tribution of  the  two  soft  limiters.  Analogously  to  Eq.  (A-68)  we  obtain 
in  this  case: 


T\  v (t) 

nXn 


/f  fer) F (sr)  r(t)  cos[v  °H 


4/  \ 5 


(A- 73) 


2 2 '> 
= •“  I y“ 
5 2 2 


( A-71) 


Here  z is  the  total  power  in  all  the  input  signals,  while  y^  and  y ^ ar< 
the  slope  parameters  of  the  two  soft  limiters. 


d.  Summary 

Using  a time-domain  analysis  approach,  general  expressions  for 
the  amplitude,  frequency,  and  phase  of  the  signals  and  intermodulation 
products  at  the  output  of  the  bandpass  limiter  have  been  obtained  for  the 
case  when  the  limiter  input  consists  of  n phase-modulated  carriers  of 
arbitrary  amplitudes  and  noise.  The  amplitudes  of  the  signals  and  cross 
products  can  be  determined  by  numerical  evaluation  of  an  infinite  integral 
for  the  different  limiter  characteristics  considered  in  Section  2.  Equa- 
tions have  also  been  obtained  to  determine  the  frequency  and  phase  of  the 
limiter  output  components.  It  has  also  been  possible  to  obtain  the  dis- 
tribution of  the  up-link  noise  at  the  bandpass  limiter  output.  Since, 
in  FDMA,  the  total  signal  power  at  the  repeater  input  is  much  larger  than 
the  noise  power,  the  distribution  of  the  noise  at  the  repeater  output  may 
be  considered  Gaussian. 

4.  Conclusions 

Analytical  expressions  for  the  nonlinear  transfer  characteristic  of 
a quasi-linear,  a soft,  and  a hard  limiter  can  be  obtained  from  an  in- 
finite integral  representation  that  relates  the  limiter  output  to  its 
input.  In  the  case  of  a soft  limiter,  the  integral  expression  is  iden- 
tical to  an  error-function  representation  commonly  used  to  characterize 
a soft  limiter.  Generally,  the  power  output  of  practical  limiters  in- 
creases linearly  (unity  slope)  at  low  input  power  levels  and  reaches 
saturation  at  large  power  levels.  Representation  of  the  limiter  transfer 
characteristics  as  a sum  of  two  soft  limiters  each  represented  by  an 
error  function  ensures  unity  slope  at  low  signal  levels,  with  an  addi- 
tional parameter  available  to  characterize  the  softness  with  which  the 
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limiter  reaches  saturation. 


The  effect  of  ideal  bandpass  limiting  on  an  arbitrary  signal  lying 
in  narrow-band,  stationary,  zero-mean  Gaussian  noise  has  been  analyzed, 
and  general  analytic  expressions  for  the  limiter  output  components  have 
been  derived  by  using  a time-domain  approach.  A major  advantage  of  the 
time-domain  method  is  that  it  preserves  the  phases  of  the  signal  and  inter- 
modulation products,  which  are  destroyed  in  the  characteristic  function 
method  generally  used  in  limiter  studies.  Expressions  for  the  desired 
signal  and  intermodulation  product  amplitudes  have  been  obtained  for  the 
case  when  the  input  consists  of  n phase-modulated  signals  of  arbitrary 
amplitudes  and  noise.  The  amplitudes  of  the  signals  and  intermodulation 
components  can  be  determined  by  numerically  evaluating  an  infinite  integral 
for  the  different  limiter  characteristics.  It  has  also  been  possible  to 
obtain  a time-domain  representation  for  the  limiter  output  noise,  which 
is  particularly  useful  for  the  evaluation  of  the  system  following  the 
bandpass  limiter.  If  the  total  signal  power  at  the  limiter  input  is  much 
greater  than  the  repeater  noise  power,  the  distribution  of  the  noise  at 
the  repeater  output  can  be  considered  Gaussian. 
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ERROR- RATE  I N VEST I GA T ION 


1.  General 

The  objective  ol  this  appendix  is  to  provide  a comprehensive  analysis 
of  the  problem  of  detection  of  n quadr iphase-modulated,  constant-envelope 
PSK/FDMA  carriers  after  transmission  through  a wideband  limiting  satellite 
transponder.  Specifically,  the  goal  of  the  study  is  to  develop  the  ex- 
pression for  the  bit-error  rate  at  the  output  of  an  KDMA  channel,  taking 
into  consideration  the  presence  of  (1)  other  KDMA  carriers  causing 
adjacent-channel  interference,  (2)  cross  products  generated  as  a result 
of  the  mixing  of  the  KDMA  signals  in  the  limiting  repeater,  and  (3)  re- 
transmitted up-link  (satellite  repeater)  as  well  as  down-link  (receiver) 
noise. 


i 


2.  Quadriphase  Detection 

A digital  quadr iphase-modulated  carrier  can  be  expressed  as: 


s(  t ) 


A cos 


tv 


0 (t) 
1 


(Ji-1) 


where  u).  is  the  carrier  radian  frequency,  and  0.  is  a constant  phase 
angle.  The  phase  modulation  c (t)  at  any  instant  can  take  one  of  lour 
values:  tt/4,  3n/4,  5tt/4,  or  7tt/4.  Dy  expanding  the  cosine,  it  can 

easily  be  seen  that  the  quadriphase  signal  can  also  be  expressed  as  the 
sum  of  two  orthogonal  biphase-modulated  signals  (each  having  half  the 
data  rate  of  the  quadriphase  modulation), 
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FIGURE  B-2  CORRELATION  RECEIVER  FOR  DETECTION  OF  A QUADRIPHASE  SIGNAL 


v.(t).  No  interference  occurs  between  the  two  channels,  because  the  two 
sinusoidal  carriers  are  in  quadrature  and  therefore  are  orthogonal. 

3.  Approach  for  Determining  the  Error  Kate 

Kor  determining  the  error  rate  of  any  I'DMA  channel,  it  will  be 
assumed  that  all  the  transmitters  in  the  system  employ  quadriphase  modu- 
lation of  the  PDMA  carriers  and  that  the  corresponding  receivers  use  the 
detection  scheme  shown  in  Figure  IJ-2.  In  practice,  a wideband  RF  band- 
pass filter  will  usually  precede  the  detector  in  Figure  B-2.  The  purpose 
of  this  filter  is  to  ease  the  problem  of  initial  receiver  synchronization 
by  removing  out-of-band  signals  and  intermodulation  products  without 
appreciably  distorting  the  desired  signal.  For  the  evaluation  of  the 
error  rate,  we  may  assume  that  receiver  synchronization  has  been  achieved 
and  consequently  the  RF  filter  is  not  required.  The  suppression  of  the 
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out-of-band  signals  and  intermodulation  products  is  accomplished  by  the 
low-pass  ( integrate-and-dump)  filters  of  Figure  B-2. 

Since  a quadr iphase-modulated  signal  can  be  expressed  as  two  orthogo- 
nal biphase-modulated  signals,  the  bit-error  rate  of  a quadriphase  system 
can  be  defined  as  the  average  of  the  bit-error  rates  of  tiie  two  biphase 
channels  in  Figure  B-2.  Because  of  symmetry,  the  bit-error  rate  of  the 
two  channels  will  be  the  same.  The  error  rate  of  a quadriphase  system 
is,  therefore,  identical  with  that  of  a biphase  system,  and  it  suffices 
to  consider  the  probability  of  error  of  either  the  upper  or  lower  channel 
in  Figure  B-2  only.  For  the  error  analysis,  the  FDMA  receivers  can  be 
represented  simply  as  single  biphase  detectors,  as  shown  in  Figure  B-3. 

BEAD  OUT 
AT 


REFERENCE 

SA  -1975-39 

FIGURE  B-3  RECEIVER  MODEL  FOR  ERROR-RATE  CALCULATION 

At  the  receiver  output  in  Figure  B-3,  the  decision  whether  a mark 
or  a space  has  been  transmitted  is  made  on  the  basis  of  whether  the  fil- 
ter output  at  the  sampling  instant  is  positive  or  negative.  If  it  is 
assumed  that  the  transmission  of  a mark  or  a space  signal  is  equally 
probable  and  that  the  noise  at  the  filter  output  has  a zero-mean  value, 
then  the  probability  of  error,  P , in  a decision  is  equal  to  the  error 
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will  be  negative  at  the  sampling  instant.  Thus,  P is  given  by 


p(z)  tlz 


(B-3) 


where  p(z)  is  the  probability  density  function  of  the  filter  output.  To 
calculate  P , one  needs  to  determine  p(z),  which  may  be  obtained  by  first 
calculating  the  characteristic  function  of  the  filter  output  and  then 
taking  the  Fourier  transform: 


p(z)  = — 
2n 


/-x  - i z 

C (?)  e ~ d? 
z 


(B-4) 


where  C (?)  is  the  characteristic  function  of  the  random  variable  z at 
z 

the  output. 

In  nonlinear  systems,  such  as  in  the  case  of  FDMA  where  the  signals 

and  up-link  noise  have  passed  before  detection  through  a limiter  having 

a nonlinear  characteristic,  this  transform  poses  a formidable  problem, 

since  the  expression  for  the  characteristic  function  is  usually  extremely 

complex.  Alternatively,  it  is  much  more  convenient  to  calculate  the 

error  probability  directly  from  the  characteristic  function.  It  can  be 

shown  that,  if  z is  any  random  variable  of  probability  density  function 

p(z)  and  characteristic  function  C (?),  its  cumulative  distribution  lunc- 

z 

tion,  P(a),  is  given  by:' 


P(a)  = p(z)  dz  = — - — I fc  (.%) 

2 -r  mLz 


■ i ? a | dc 


( B-5 ) 


probability  in  either  a mark  or  a space  signal.  The  probability  of  error 
in  detecting  a mark  is  equal  to  the  probability  that  the  filter  output 
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where  I denotes  that  only  the  imaginary 
m 

par  t 

has  to  be  taken. 

Compari son 

- j 

of  Eqs.  (B-3)  and  (B-5)  shows  that 

CO 

1 1 

f 1 

- 

: 

P = P(0)  = 

I 

c (-: 

(B-6) 
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Thus,  the  approach  for  calculating  P for  any  FDMA  channel  will  he 

to  determine  C (£)  at  the  channel  output  for  a mark  transmission  and  then 
Z 

to  perform  the  integration  in  Eq.  (B-6).  It  should  be  noted  that  Eq.  (B-6) 

is  a general  result  that  can  be  used  to  determine  P for  a linear  as  well 

e 

as  a nonlinear  channel  of  arbitrary  transfer  characteristic. 

4.  Determination  of  C7(-) 

To  derive  the  mathematical  expression  for  C (£)  at  the  receiver  out- 

z " 

put  in  Figure  B-3,  it  will  be  necessary  first  to  obtain  the  expression 
for  the  receiver  output  z at  the  sampling  instant.  To  achieve  this  we 
write  an  expression  for  the  signals  and  noise  components  present  at  the 
receiver  input  and  then  develop  the  expression  for  the  output. 

Appendix  A shows  that  when  the  input  to  the  satellite  repeater  con- 
sists of  n constant-envelope  quadr iphase-modulated  FDMA  carriers  and  up- 
link noise,  the  repeater  output  will  contain,  in  addition  to  the  signals 
and  noise  components,  intermodulation  products  resulting  from  the  mixing 
of  the  signals  in  the  limiter.  All  these  components  will  also  be  present 
at  each  receiver  input,  but  their  amplitudes  will  be  modified  because  of 
the  power  loss  on  the  down-link  transmission  and  the  amplification  in  the 
receiver  front  end.  In  addition  to  these  signals,  received  from  the 
satellite,  down-link  noise  will  be  present  at  the  input  to  the  receiver. 

th 

Thus,  in  an  FDMA  system  the  input  to  the  i receiver  will  comprise 
the  components  considered  below. 
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The  Desired  Signal  Component 

s (t)  = a [cos  u)  t f (t)  9 "1 

1 i L i i i-J 

a 

= ~=r  [u  (t)  cos (x  t h 0 ) - v (t)  sink  t 
2t_i  \i  i / i Vi 
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(13-7) 


where  u (t)  and  v (t)  are  independent  binary  (1)  data  modulation  of  the 
i i 

orthogonal  carriers,  and  Q is  a random  phase  angle. 
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The  Undesired  Signal  Components 
n 


s ( t)  = 
k 


a cos  !'x  t 0 ( t ) 
k L k k 


• 


(B-8) 


k = l 
k^  i 


where  ^ ( t ) represents  quadriphase  modulation  of  the  signals,  and  9 is 
the  random  phase  of  the  carriers. 


The  Intermodulation  Components 
m 

I(t)  = ^ 

7=T 


hi 

> b cos  |oi  t 1-  ■■  ( t)  ^ , 

! L .1  .1  .1  J 


( B-9 ) 


where  ■;  (t)  represents  the  phase  modulation  of  the  intermodulation 

.1 

products;  n.  is  the  random  phase  angle;  and  m is  the  number  of  cross 
products  falling  in  the  repeater  bandwidth. 

Expressions  for  determining  the  amplitudes,  phases,  and  frequencies 
of  the  signals  and  intermodulation  products  at  the  output  of  the  satellite 
repeater  are  given  in  Appendix  A.  The  expressions  for  the  signal  ampli- 
tude and  the  intermodulation  products  at  the  receiver  input  can  be  ob- 
tained from  their  corresponding  expressions  at  the  satellite  repeater 
output  by  taking  into  consideration  the  losses  on  the  down  1 ink  and  the 
gain  in  the  receiver.  The  frequencies  and  phases  of  the  intermodulation 
products  in  Eq.  (B-9)  are  dependent,  respectively,  on  the  frequencies 
and  phases  of  the  signals.  The  number  of  intermodulation  products  falling 
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in  the  repeater  bandwidth  m is  solely  a function  of  the  number  of  IDMA 


signals  at  the  repeater  input  ami  their  frequencies.  lor  the  calculation 
of  the  error  rate  we  assume  a frequency  assignment  for  the  signals,  so 
the  number  of  in-hand  intermodulation  products  m,  their  frequencies  x., 
and  their  phases  •![  ( t)  , a can  be  determined  from  the  knowledge  of  the 

.i  i 

frequencies  and  phases  of  the  signals  producing  these  cross  products. 

The  Noise  Components 

The  receiver  input  noise  consists  of  the  sum  of  the  retransmitted 
up- link  noise  and  the  Gaussian  down- link  or  receiver  front-end  noise. 
Appendix  A shows  that  the  distribution  of  the  up-link  noise  at  the  re- 
peater output  will  be  to  a good  approximation  Gaussian,  since  the  total 
FDMA  signal  power  at  the  repeater  input  is  much  larger  than  the  up-link 
noise  power.  Thus  the  distribution  of  the  retransmitted  up-link  noise 
at  the  receiver  input  will  also  be  Gaussian;  only  its  variance  relative 
to  that  at  the  repeater  output  will  be  different.  This  can  be  determined 
by  taking  into  account  the  down-link  losses  and  the  amplification  in  the 
receiver.  Since  both  the  up-link  and  the  down-link  noise  at  the  receiver 
input  are  independent  and  zero-mean  Gaussian,  the  composite  noise  will 
also  be  Gaussian  with  zero  mean  and  variance  equal  to  the  sum  of  the 
variances. 

The  input  to  the  i**1  receiver  in  Figure  13-3  can  be  expressed  as 
a 

r(t)  = [u  (t)  cos|x^t  + - v^Ct)  sin^x^t  F 0 . j]  -F  n(t) 

1 \ a cos  ^0)  t t-c  (t)  P 0 1 
/ ^ k - k k k-l 

k 

P £ b cos  [u)tP’i.(t)  O'”!  . ( B- 10) 
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Here,  n(t)  denotes  the  composite  receiver  input  Gaussian  noise  (sum  of 
retransmitted  up-link  and  down-link  noise). 


The  receiver  output  z in  Figure  13-3  is  obtained  by  cross  correlating 

the  input  r(t)  with  the  synchronized  receiver  reference  2 cos| x t Q ] 

i i 

over  the  bit  interval  T of  the  desired  signal, 
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(Note  T is  equal  to  2 R where  R is  the  bit  rate  of  the  quadriphase 
i i i 

signal.)  Substituting  Kq.  (13-10)  into  Kq.  (13-11)  and  assuming  that  a 

mark  was  transmitted  [u  (t)  = 1|,  the  output  is  given  by: 
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The  carrier  with  modulation  v (t)  in  Kq.  (13-10)  will  not  produce  any  out- 

i 

put,  since  it  is  orthogonal  to  the  receiver  reference.  The  parameter 
in  Eq.  (13-12)  denotes  the  receiver  output  noise.  Since  the  receiver  input 
noise  n(t)  is  assumed  to  be  Gaussian  and  the  receiving  system  is  linear, 
the  distribution  of  the  output  noise  will  also  be*  Gaussian  with  zero 
mean  and  variance  given  by 


' / 


G ( f ) df 


( 13-1:3) 


.» 1 


where  G(f)  is  the  spectral  density  of  the  receiver  input  noise,  and 
H(f)  is  the  transfer  function  of  the  integrate-and-dump  filter  in  figure 

B-3. 


By  expanding  the  cosines  in  Kq.  (B-12)  and  rearranging  the  output, 
z can  be  written  as: 


z = — - - f T1 

r?>  ' 


A cos  0 
ik  ik 


EB  cos  a 

i.i  i.J 


(B-14) 


where  the  amplitudes  and  phases  of  the  undesired  signals  and  intermodula- 
tion products  at  the  receiver  output  are  defined  as: 


A = a 
ik  k 


B = b 
i.j  J 


y ik  ik 


2 2 

3.  . fV.. 
i.l  iJ 


(B-15) 


) =0  - 0.  t arc  tan  y /f3 

ik  k l ik  ik 


r v = a - 9.  J arc tan  y.  . 7 3 . . 

ij  .1  i ij  iJ 


(B-16) 


The  quantities  3 and  y in  Eqs.  (B-15)  and  (B-16)  represent  the  elfect  ot 
filtering  on  the  signals  and  cross  products.  They  are  given  by: 
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ik 


= Ti  J 
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• i 

cos[| 
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X - X V 

ik  i/ 

i a,  ( t )1 

k J 

dt 

(B-17a) 
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The  S and  V 

ik  ik 


as  we  1 1 as 


ij 


and  y. 


i.J 


are  random  variables,  be- 


cause they  are  dependent  on  the  random  modulation  C (t)  of  the  FDMA 

k 


It  should  be  remembered  that  the  phase  modulation  Oj  (t)  of 

j 


carriers. 

the  cross  products  is  a function  of  the  modulation  e(t)  of  the  signals. 


Equation  (B-14)  represents  the  output  of  the  receiver  when  its  input 
consists  of  signals,  cross  products,  and  noise.  To  calculate  the  error 
probability  by  Eq.  (B-6),  we  need  the  characteristic  function  C (?)  of 

7.  ~ 

the  receiver  output  z.  The  first  two  terms  in  Eq.  (B-14)  represent  the 
sum  of  the  desired  signal  and  the  contribution  of  the  receiver  input  noise 
at  the  output.  The  last  two  terms  containing  the  summations  represent 
the  contributions  of  the  undesired  signals  and  the  intermodulation  products. 
Since  the  first  two  terms  are  statistically  independent  of  the  last  two, 

C (§)  can  be  expressed  as  the  product  of  the  characteristic  function  of 
the  signal  plus  noise  and  the  characteristic  function  of  the  sum  of  the 
undesired  signals  and  the  cross  products: 

•>  •> 

i?a.  2 - '75  2 

C {")  = e 1 • e ‘ • f < - ) . (B-1H) 


Here,  C(|)  denotes  the  characteristic  function  of  the  sum  of  the  unde- 
sired signals  and  the  intermodulation  components  in  Eq.  (B-14). 

The  angles  Q and  a in  Eq.  (B-14)  are  uniformly  distributed  and 
ik  i.j 

are  also  statistically  independent.  The  amplitude  A ^ of  any  one  unde- 
sired signal  at  the  i^  receiver  output  is  dependent  only  on  the  random 

modulation  ® (t)  of  that  signal.  But  the  amplitude  B ol  any  given 
k i.j 


5(1 


cross  product  is  dependent  on  the  phase  modulations  of  the  signals  that 

generate  that  cross  product.  Averaging  first  over  the  distributions  of 

0 and  ry  gives: 
ik  i.j 


C(£)  = E 


TT  • ?)  • TTJo(uu  ■ !) 


(B-19) 


where  EL  I denotes  the  expectation  of  A 


ik 


random  phase  modulation  0 of  the  signals. 

k 

into  Eq.  (B-18)  yields: 


and  B with  respect  to  the 
i.j 

Substitution  of  Eq.  (B-19) 


• n? 


C (-)  = e 

7. 
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( B-20 ) 


In  the  special  case  when  the  EDMA  signals  are  unmodulated  I 0 ( t ) 0|, 

k 

all  A and  B will  be  constants  and  we  would  not  need  to  evaluate  the 

ij  i j 

expectation  in  Eq.  (B-20). 


In  the  derivation  of  Eq.  (B-20),  it  is  implicitly  assumed  that  the 

amplitudes  of  the  signals  and  cross  products  (a.,  a , b.)  at  the  receiver 

i k ,j 

input  are  constants.  Strictly,  because  of  the  presence  of  up-link  noise 
at  the  limiter  input,  these  amplitudes  will  be  random  variables  dependent 
only  on  the  amplitude  distribution  (Rayleigh  distribution)  of  the  up-link 
noise.  Thus,  Eq.  (B-20)  should  also  be  averaged  over  the  amplitude  dis- 
tribution of  the  up-link  noise.  Since,  in  EDMA,  the  total  signal  power 
at  the  limiter  input  is  much  larger  than  the  input  noise  power,  the  fluc- 
tuation in  the  amplitudes  of  signals  and  cross  products  caused  by  the 
presence  of  noise  will  be  small  and  can  be  neglected  without  producing 
significant  error  in  the  determination  of  C (').  Thus,  instead  of 

averaging  Eq.  (B-20)  we  assume  that  the  amplitudes  a , a,  , and  a are 

1 k ,j 

the  expected  values,  (averaged  over  the  amplitude  distribution  ol  the 
input  noise).  They  are  then  identical  with  the  h parann  ters  discussed 


in  Appendix  A,  with  the  parameter  cy  modified  to  take  into  account  the 


gains  and  losses  from  the  limiter  output  to  the  receiver  input.  This 
approximation  includes  in  the  analysis  the  suppression  and  power-sharing 
effects  of  the  up-link  noise  and  neglects  the  slight  fluctuation  caused 
by  noise  in  the  amplitudes  of  signals  and  cross  products  at  the  limiter 
output. 


5 . Determination  of  P,, 

Substitution  of  the  imaginary  part  of  Eq.  (li-20)  in  Eq.  (B-6)  yields 

the  following  expression  for  the  bit-error  probability  at  the  output  of 
th 

the  i channel: 


( B-21 ) 


It  can  be  seen  that  in  Eq.  (B-21)  the  sin  E E term  is  a maximum  at  E = 0 
and  decreases  rapidly  as  § increases.  The  same  is  true  for  the  exponen- 
tial term.  The  value  of  J ( ) is  unity  for  E = 0 and  less  than  unity  for 

o 

all  other  values  of  As  the  number  of  terms  in  the  product  increases, 

the  value  of  the  product  becomes  small,  except  when  all  terms  are  near 
unity,  i.e.,  when  § is  near  zero.  Since  the  integrand  of  Eq.  (B-21)  is 
the  product  of  the  three  terms,  which  decrease  with  increasing  §,  the 
major  contribution  to  the  integral  will  stem  from  the  values  of  these 
terms  near  § 0.  In  the  neighborhood  of  F = 0,  we  can  expand  the  Bessel 

functions  in  a power  series  and  use  the  following  approximations: 
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TT  j (a 
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TT  • 5) 
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= 1 - •*—  E 
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Aik  ^ 

/ 1 o \ i k 
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Z-f  2 ' A 
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I terms  in  § and  higher 


where 


(B-22) 


> 

( li-23  ) 
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(D-24) 


are  respectively  the  average  signal  and  cross-product  powers  at  the  re- 
ceiver output.  The  summations  in  Eq.  (B-22)  represent  the  total  power 

contribution  of  all  the  undesired  signals  and  the  intermodulation  products 

th 

at  the  output  of  the  i receiver. 

The  left-hand  side  of  Eq.  (B-22)  represents  the  product  of  the  charac- 
teristic functions  of  independent  sinusoids.  If  the  number  of  sinusoids 
is  large,  the  resultant  characteristic  function  can  be  approximated  by 
the  characteristic  function  of  Gaussian  noise  |Eq.  (B-23)|,  which  is 
.justified  by  the  central  limit  theorem.  lor  as  few  as  four  signals,  we 
have  in  Eq.  (B-22)  a product  of  27  Bessel  functions  (three  signals,  24 
third-order  cross  products).  Eor  ten  signals,  the  number  of  Bessel  func- 
tions increases  to  459.  So  it  can  be  seen  that  the  approximation  of  Eq. 
(B-23)  will  improve  as  the  number  of  signals  increases. 


Using  the  approximation  of  Eq,  (15-22)  in  Eq.  (B-21)  and  evaluating 


the  integral  yields: 


p.»  ■ 5 1 ■ err(»//T)]  3 E-n.  E'Vi 


- D . 


(15-25) 


2 2 2 

where  p = a '2'  . On  the  other  hand,  substituting  Eq.  (15-23)  into  Eq. 
i i 

(B-21)  and  performing  the  integration  gives: 


Pei  -K‘  - 'r,(ri/AH 


( B-26 ) 


whe  re 


2 

r = 
i 


2 2 
a /2: 
i i 


Vp  2 Vp  2 

/ j ik  i / j l.j  i 


1 t 


( 15-27 ) 


k J th 

represents  the  equivalent  SNR  at  the  output  of  the  i channel. 


It  can  be  seen  from  the  numerical  tables  that  lor  small  x the  value 

of  J (x)  lies  between 
o 


2 2,  , 

x -x  4 

1 - — < J ( x)  <.  e 
4 o 


Thus,  the  approximation  of  Eq.  (B-22)  will  yield  a lower  bound  on  the 

error  probability  (Eq.  ( B—  25)1,  while  the  exponential  approximation  will 

give  the  upper  bound  |Eq.  (B-26)|.  The  true  probability  of  error  will 

lie  between  these  two  values.  As  the  number  of  signals  n increases,  the 

2 

number  of  third-order  cross  products  increases  as  n (n  - 1)  2,  and  the 

total  number  of  Bessel  functions  involved  in  Eq.  (15-21)  increases  as 
2 

(n  - 1)  (-  n (n  - l)  2.  Thus,  from  the  central  limit  theorem,  it  is 

reasonable  to  expect  the  true  error  rate  to  be  closer  to  the  upper  bound 

given  by  Eq.  (B-26).  We  therefore  use  this  upper  bound  to  determine  the 

th 

error  rate  at  the  output  of  the  i EDMA  channel. 


The  output  SNR  of  Eq.  (B-27)  can  also  be  expressed  as 


where 


and 


(11-28) 
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p = a 2 = SNR  lor  the  desired  signal 
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2 2,2  th 

p = a 2 = SNR  for  the  k si  anal 
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o = b /2c  = mtermodulat ion-to-noise  ratio  for  the  i cross 
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s = E^  . 1 
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( B- 29 ) 


The  parameters  S and  S represent  the  effect  of  filtering  in  the  i 
ik  ij  ^ 

receiver  on  the  k interfering  link  and  on  the  j cross  product, 
respectively. 


tit 


The  approximation  of  Kq.  (B-23)  amounts  to  modeling  the  composite 
distribution  of  the  sum  of  all  the  undesired  EDMA  signals  and  the  cross 
products  as  equivalent  Gaussian  noise  at  the  output  of  the  receiver. 

This  approximation  is  not  the  same  as  assuming  their  distribution  to  be 
Gaussian  at  the  receiver  input.  If  care  is  not  exercised,  the  latter 
approximation  will  lead  to  errors  in  the  calculation  of  the  variance  of 
the  Gaussian  distribution  at  the  receiver  output.  Kor  example,  consider 
the  contributions  of  the  two  third-order  cross  products, 


58 


and 


A cos  (x 

X,  - x )t 

+■  rv 

tv  " tv , 

L\  i 

l 2/ 

i 

1 2-J 

A f-os[(:.  I -a>2  - U)jl  • rv.  ' rv2  - nrj 

to  the  receiver  output.  These  are  generated  as  a result  oi'  the  mixing 

of  the  desired  signal  at  the  frequency  f with  two  other  signals  at  fre- 

t 

quencies  f and  1^.  If  these  cross  products  are  considered  independent 
at  the  receiver  input  (this  is  required  by  the  Gaussian  assumption), 
their  contributions  to  the  receiver  output  can  be  obtained  by  summing  the 
power  contribution  of  each,  taking  into  account  their  power  spectra  and 
the  transfer  function  of  the  integrate—  and-dump  filter. 

However,  more  accurately,  the  contribution  of  these  two  cross  products 
is  obtained  by  cross  correlation  with  the  receiver  reference.  Multiplica- 
tion of  both  the  cross  products  with  the  receiver  reference  will  yield  a 

single  baseband  waveform,  2A  cos|(  j-  - x ) t rr  - n-  I , of  twice  the 

12  12’ 

amplitude,  because  the  two  cross  products  add  coherently  in  amplitude. 

The  power  contribution  at  the  receiver  output  can  then  be  obtained  by 
integrating  the  power  spectrum  multiplied  by  the  square  of  the  magnitude 
of  the  transfer  function  of  the  filter.  This  more  accurate  procedure 
leads  to  a 3-dB-larger  value  for  the  output  power  contribution  from  these 
two  cross  products  than  the  one  obtained  by  considering  them  as  inde- 
pendent sinusoids  at  the  receiver  input. 

The  time-domain  approach  takes  into  account  any  coherency  between 
cross  products  at  every  stage  of  the  receiver;  it  therefore  provides  an 
accurate  determination  of  their  contribution  to  the  receiver  output. 

The  approximation  of  Eq.  (ii-23)  is  based  on  determining  the  contribution 
of  all  the  undesired  FOMA  signals  and  cross  products  at  the  receiver  out- 
put, taking  into  account  their  coherency,  and  then  approximating  them  as 


5i> 


Gaussian  noise. 


Determination  of  and  Sj 


Consider  first  the  determination  of  S : 

i k 


s = K/32  t-  v2  ) , 

ik  \ ik  ik* 


^ik  = 


f f cosK 


0)1  ' 0 ( t)  (It 

lv  k — * 


r f sin[Av  1 Vl)] 


( B-30 ) 
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Here,  the  expectation  has  to  be  evaluated  over  the  distribution  of  © (t) 

k 

during  the  bit  interval  0 < t <■  T . 

i 

2 

'•rom  Kq.  (B-30)  it  can  be  seen  that  the  expectation  of  g ^ will  be 

the  same  as  that  of  v . Thus,  we  can  write 

ik  ’ 


S = 2 • K E 

ik  V 


(B-31) 


'I ne  expectation  of  g.  can  be  calculated  as  follows: 

1 K 


T T 

• i r i 


E(?iJ = ^ e!cos[aVi  1 °k(ii)]  • cos[%t2  * \(i2)]!  ciiidt2 


l o o 


T T 

• i r i 


E , cos 


■K(s  - s)  \(iJ-  %(•■)] ! "vs 


l o o 

T T 

r i r i 


Ro  ( 1 2 - ll)  COSK(l2  " tl)]dtl‘,t2  ’ (B-32) 
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where  R0  (t  - 1 ) represents  the  autocorrelation  function  of  the  quadri- 
k 2 1 

phase  modulation  c (t).  The  Last  line  of  Eq.  (15-32)  is  obtained  by  ex- 
it 

panding  the  cosine  in  the  previous  line  and  using  the  relationships 

uc  (i2  - s)  = E[sin  \(ii)  • sin  \(t2)] = ECcos  °k(ii)  • cos  \(i2)] 


E[sin  • JtJ  • cos  *k(t2)]  = K[sin  0R(tJ 


0 . (15-33) 


Defining  the  impulse  response  of  the  integrate  and  dump  filter  in 


figure  15-3  as: 


h ( t)  = 1 'T  for  0 t • T 

i i i 


= 0 t > T. 


we  can  express  Eq.  (15-32)  as: 


00  00 

E(3ik)  = J | \([2  - s)  COS  AXu(l2  “ ll)  • hi(‘l)  • hi(s)  “V 


(15-31) 
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Replacing  t by  t^  i t in  the  above  expression  yields: 

00  00 

E(3^)  = [ I \{T)  • cos  AV  • hi(ti)  • hi(ti  ' T)  dt! 

O —oo 

If  the  filter  autocorrelation  function  is  defined  as 


dr  . (15-35) 


R (t)  = 
h 
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(B-36) 


we  obtain,  finally,  the  relationship: 


PJ  - f % 


(t)  • R (t)  • cos  Au)  T cii 
h k 

k i 


Substituting  Eq.  (15-37)  in  Eq.  (15-31),  we  obtain: 


(15-37) 


- « f ■ 


(r)  • R (t)  ■ cos  Ait  T dT 
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(15-38) 


It  can  be  shown  that  the  filter  autocorrelation  Rh  (t)  is  a triangle 

"i 

given  by: 


, <t>  * f (‘  - ¥) 


for  T - T 


(B-39) 


The  autocorrelation  function  of  the  quadriphase  signal  modulation  Ko.(t) 

K 

can  be  determined  as  follows: 


(t)  = E (cos  © (t)  • cos  0 (t  r )1 
1-  k k -I 


= — E [u  (t)  • u (t  t)1 
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( 15-  ’0  ) 
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where  u (t)  is  a binary  ( 1)  random  waveform  whose  autocorrelation  tunc- 

k 

tion  is  represented  by  Ru  (t).  Since  the  autocorrelation  function  of  a 

k 

binary  random  waveform  is  a triangle,  It  (T)  can  be  expressed  as: 

uk 


R (r) 
n 
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■(■-?)  ■ 


for  It  I - T 
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(B-41) 


where  1/T  is  the  bit  rate  of  the  biphase  modulation.  The  parameter  S 

k i k 

can  now  be  evaluated  with  the  aid  of  Eqs.  (B-39)  through  (B-41). 

Similarly,  for  S . we  obtain  the  following  expression  analogous  to 
Eq.  (B-38) : 


S =2  I It  (t)  • It  (t)  cos  /\t  t dr 
BI  h.  j 
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(B-42) 


Ax . = x . - x 
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and  R...  . (t)  is  the  autocorrelation  function  of  the  phase  modulation  of  the 


intermodulat ion  produc t 


R (t)  = E cos  (t)  • cos  ' (t 

L J M 


( B-43 ) 


The  evaluation  of  FU  (t)  is  best  described  by  an  example.  Suppose  we 
want  to  determine  the  au  ^correlation  function  of  a third-order  cross 
product  having  the  phase  modulation  't(t)  = © (t)  f ©,^(t)  - ®^(t).  This 
cross  product  is  generated  as  a result  of  mixing  three  input  signals  in 
the  limiter  with  the  phase  modulations  ®^(t),  ©,,(t),  and  ©^(t),  respec- 
tively. In  terms  of  signal  modulations,  Eq.  (B-43)  can  be  expressed  as: 


(15-44) 
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R (t)  = e|cos[©  (t)  5 0 (t)  - <s>  (t)~l 
t . | L i 2 4 J 


cos  ( t t t)  c^(t  t t)  - ®.{(t  * T)j 


By  expanding  the  cosines  and  using  the  statistical  independence  between 
the  phase  modulations  of  the  signals,  it  can  be  shown  that 


R (r)  = 4R  (t)  • R (t)  • R (t) 
© e © 

j 1 2 3 


(IS-45) 


where  the  R (t)'s  are  the  autocorrelation  functions  of  the  random  quadri- 
o 

phase  modulation  of  the  signals.  Each  R (t)  can  be  expressed  in  terms  of 

© 

a biphase  autocorrelation  function  R (r),  as  in  Eq.  (B-38).  He  thus 

u 

obtain 


R (t)  = ~ R (t)  • R (t)  • R (t) 
•'r  2 u u u 

j 12  3 

where  each  R (t)  is  a triangle  defined  as: 


( 15-46) 


R (t) 
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-i-M  , 

rp 


for 


i = 1,2,3  . ( 15-47 ) 


Equation  (K-42)  can  now  be  evaluated  with  the  aid  of  Eqs.  (B-46)  and  (B-47), 

Using  this  approach,  we  can  determine  the  power  contribution  of  any 
desired  intermodulation  product  at  the  receiver  output.  A special  case, 
however,  occurs  in  dealing  with  the  cross  products  generated  as  a result 
of  the  mixing  of  the  desired  signal  with  other  signals.  For  example, 
consider  the  case  where  we  want  to  determine  the  power  contribution  of  a 
third-order  intermodulation  product  that  has  the  modulation 


( t)  = « ( t)  h © ( t)  - 0 ( t) 
1 2 j 
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1 


Since  the 


where  C (t)  is  the  phase  modulation  of  the  desired  signal. 

receiver  will  be  synchronized  with  the  desired  signal,  i.e.,  with  C (t), 

the  output  power  for  this  cross  product  will  be  independent  ol  o (t)  and 

i 

will  depend  only  on  the  random  switching  of  c^(t)  and  c^(t)  during  the 

integration  interval  of  the  matched  low-pass  filter.  The  autocorrelation 

function  R,  (t)  in  this  case  will  be  determined  by  the  autocorrelation 

functions  of  o and  © : 

1 2 


R,  (t) 
' j 


E j cos  t ( t ) ©_  ( 

I — 1 2 


2 R (t)  • R (t) 
© 0 
1 2 


= - R (r)  • R (t) 


■ cos  © ( t 

t) 

0(1  T)1 

-1  L 1 

2 -1 

(13-48) 


Similar  expressions  can  be  derived  for  treating  cross  products  of  other 
orders  that  are  generated  by  mixing  of  the  desired  signal  with  other  sig- 
nals in  the  limiter. 


7 . Summary 

By  using  a time-domain  analysis  approach,  a general  expression  [Eq. 
(B-21)|  has  been  obtained  for  determining  the  bit-error  probability  in 
detecting  a desired  quadr iphasc-modulated  RDM A signal  at  the  receiver, 
when  n independent  FDMA  signals  are  transmitted  simultaneously  through  i 
limiting  satellite  transponder.  The  error-rate  expression  has  been  de- 
rived  under  the  assumption  that  small  fluctuations  in  the  amplitudes  ot 
the  signals  and  cross  products  caused  by  the  random  envelope  ol  thi 
limiter  input  noise  can  be  ignored  and  that  the  distribution  oi  the  up- 
link noise  at  the  output  of  the  repeater  can  be  considered  t.aussi  a. 

In  an  FDMA  system,  these  assumptions  are  permissible  under  the  condition 

j 

I 


3 


that  the  total  power  in  all  the  signals  at  the  repeater  input  is  much 
larger  than  the  power  in  the  up- link  noise. 


Exact  evaluation  of  Eq.  (B-21)  is  very  difficult  because  ol  the 
problem  involved  In  determining  the  expectations  of  the  rec<  iver  output 
amplitudes  of  the  interfering  KDMA  signals  ind  the  cross  products  caused 
by  random  modulation  of  the  carriers.  However,  it  has  been  possibli  to 
obtain  the  upper  and  lower  bounds  within  which  tilt  error  rate  is  expt  c ■ <t 
to  lie.  When  the  number  of  FDMA  s if nals  is  larg<  error  I ex- 

pected to  be  close  to  the  upper  bound. 
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Appendix  C 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM  SYSCO N 

1 . Gene  ra 1 

The  computer  program  SYSCON  has  been  developed  to  model  an  FDMA 
satellite  communication  system  and  to  investigate  and  optimize  its  per- 
formance in  operation  with  a mix  of  user  terminals  of  differing  charac- 
teristics, such  as  data  rates,  transmitter  ERPs,  and  receiver  sensitivi- 
ties (G/T).  The  capability  to  improve  system  performance  is  achieved 
through  optimization  of  a norm  representing  the  figure  of  merit  or  a 
measure  of  the  system's  communication  performance. 

Two  norms  have  been  selected  for  optimization: 


n 

N = 
a / j el 
i = 1 


(C-l) 


and 


n / n 

XX  ■ RyXX 

i = l ' i-i 


N = 

b / j ei 
i = 1 


(C-2 ) 


where 


error  probability  of  the  i channel 


R = data  rate  of  the  i 
i 


th 


channe 1 . 


The  norm  N treats  all  links  equally  and  does  not  take  into  consideration 
a 

the  fact  that  some  links  may  carry  information  at  higher  data  rates  than 


o t hers . 


This  fact  is  taken  into  account  in  norm  N , which  represents  the 

b 


till 
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sum  of  error  probabilities  weighted  by  the  data  rates.  Either  of  these 
two  norms  can  be  selected  for  system  optimization  in  SYSCON. 

2.  Model  of  the  Satellite  Repeater 

The  model  of  the  satellite  repeater  used  in  SYSCON  is  shown  in  Fig- 
ure C — 1 . It  consists  of  a receiving  antenna  followed  by  a bandpass  filter, 
a high-gain  preamplifier,  a limiter,  a bandpass  filter,  a TUT  amplifier, 
and  finally,  the  transmitting  antenna.  Both  the  preamplifier  and  the  TUT 
amplifier  are  assumed  to  be  linear,  so  all  the  nonlinearity  existing  in 
the  repeater  can  be  attributed  to  the  presence  of  the  bandpass  limiter 
sandwiched  between  the  two  amplifiers.  Both  the  bandpass  filters  are 
assumed  to  be  ideal.  The  input  filter  is  wide  enough  to  pass  the  signals 
with  negligible  distortion  and  limits  the  input  noise  to  a bandwidth  that 
is  small  compared  to  the  frequencies  of  the  signals.  The  bandpass  filter 
following  the  limiter  confines  the  limiter  output  spectrum  essentially 
only  to  the  fundamental  band  of  the  input  signals.  In  practice,  the  re- 
peater will  also  contain  a frequency  translator  that  ideally  has  no  effect 
on  the  system  performance.  In  SYSCON  the  up-link  frequencies  are  shifted 
down  by  an  amount  Af  (A f = 0.725  GHz)  prior  to  transmission  through  the 
limiter.  This  shift  translates  the  up-link  signal  frequencies  to  their 
correct  values  for  down-link  transmission. 

3.  Limiter  Characteristic 

The  measured  power  input/output  transfer  characteristic  of  a tunnel- 

diode  bandpass  limiter  is  shown  in  Appendix  A,  Figure  A-7.  Both  the  input 

and  the  output  power  levels  are  normalized  with  respect  to  the  maximum 

power  output  of  the  limiter  P . The  dashed  line  in  Figure  A-7  shows 

omax 

the  normalized  power  transfer  characteristic  of  a linear  amplifier  of 
unity  gain. 
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To  approximate  closely  the  transfer 


of  practical 


limiters  such  as  that  shown  in  Figure  A-7,  the  analytical  representat ion 
must  provide  unity  slope  at  low  signal  levels,  with  an  additional  param- 
eter available  to  characterize  the  softness  with  which  the  limiter  reaches 
saturation.  This  can  be  achieved  by  describing  the  limiter  voltage  char- 
acteristic by  a sum  of  two  error  functions  [Appendix  A,  Eq . (A— 32 ) ] : 


(c-:n 


where  Of  is  the  limiting  level.  The  parameters  and  \ can  be  chosen 
appropriately  to  approximate  the  measured  power  transfer  characteristic 
(such  as  in  Figure  A-7)  of  practical  limiters. 


a.  Normalization 


It  is  convenient  first  to  normalize  the  maximum  output  power  of 

the  limiter  to  0 dBW  for  the  determination  of  the  two  limiter  parameters 

V and  V (the  subscript  0 refers  to  the  values  with  respect  to  O-dBW 
10  20 

power  output).  The  values  of  v and  \ can  then  be  calculated  for  an 
arbitrary  maximum  limiter  power  output  in  terms  of  and  Y^ (J ‘ This  pro- 

cedure is  equivalent  to  incorporating  a normalized  bandpass  limiter  between 
the  preamplifier  and  the  TUT  amplifier  in  Figure  C-l.  The  resulting  model 
is  shown  in  Figure  C-2.  The  normalized  limiter  consists  of  an  amplifier 

of  gain  G preceding  the  limiter,  followed  bv  another  amplifier  of  recip- 
o 

rocal  gain  1/G.  The  gains  of  these  two  amplifiers  are  such  that  the 
power  levels  at  the  output  of  the  preamplifier  and  at  the  input  to  the 
TWT  amplifier  are  the  same  as  in  Figures  C-l  and  C-2. 


equations 


Appendix  A shows  that  Y 
[Eqs.  (A-34)  and  (A-38)l 


and 


V 

2 


can  be  determined  by  solving  two 
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(C-4  ) 


\ • iFi(i’2' 


F I-,  2, 
1 l|2 


(C-5) 


Equation  (C-4)  is  obtained  by  requiring  that  at  low  input  levels  the  output 

power  level  should  increase  linearly  with  the  input  power,  while  Eq . (C-5) 

specifies  an  output  power  P at  an  input  power  P . Since,  for  the  nor- 

ol  io 

2 2 . r 

malized  limiter  the  maximum  output  power  is  0 dBW,  we  have  a = tt  /8  lEq. 

(A-29)J.  Specifying  an  output  power  P at  an  input  power  P = P =0 

ol  io  omax 

dBW  (this  corresponds,  for  example,  to  the  normalized  output  P /P  = -4 

ol  omax 

dB  at  0-dB  normalized  input  power  level  for  the  limiter  characteristic  shown 
in  Figure  A-7),  we  can  determine  the  limiter  parameters  V and  v (refer- 
enced to  0-dBW  output)  by  solving  Eqs.  (C-4)  and  (C-5). 

In  Figure  C-2  the  limiter  is  preceded  by  two  amplifiers  of  total 

gain  G • G . If  we  now  define  P to  be  the  power  at  the  input  to  the  pre- 
a o i 

amplifier,  the  power  at  the  input  to  the  limiter  is  then  G • G P . This 

a o i 

corresponds  to  simply  dividing  the  values  of  Y, „ and  y by  G • G , as  can 

10  20  a o 

easily  be  seen  from  Eq.  (C-3) . Thus,  in  terms  of  the  power  level  at  the 
input  to  the  preamplifier  in  Eqs.  (C-l)  and  (C-2),  the  limiter  parameters 


Y and  Y are  given  by 
12 


2 10 
v = 

1 C.  • G 

a.  o 


2 C.  • G 
a o 


(C-6) 


2 2 2 2 

Here  V and  y are  the  values  of  v and  y for  G = G =0  dB. 
10  20  1 2 a o 


b.  Determination  of  G 


In  Figure  C-2  the  maximum  power  output  of  the  bandpass  limiter 
is  0 dBW,  which  is  achieved  at  large  input  power  levels.  This  power  output 
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will  yield  the  maximum  satellite  effective  radiated  power  P 


Since 


P is  given  by 
max 


ma  x 


P (dBW)  = O(dBW)  - G (dB)  + G (dB)  + G (dB) 
max  o twt  st 


we  obtain,  in  ratios  for  G , the  expression 

o 


G = 
o 


G 

st 


max 


< C — 7 ) 


where 

G = satellite  transmitting  antenna  gain 
st 

G = gain  of  the  TUT  amplifier, 
twt 

In  SVSCON,  v and  . are  determined  from  Eqs.  (C-6)  and  (C-7). 

1 2 

First,  Y and  V are  calculated  in  the  subroutine  PINOUT  by  solving 
10  20 

Eqs.  (C-4)  and  (C-5)  and  normalizing  the  maximum  power  output  of  the 

limiter  to  0 dBW.  This  corresponds,  for  example,  to  approximating  the 

power  transfer  characteristic  of  Figure  A-7  normalized  bv  a maximum  power 

output  of  0 dBW.  The  ratio  P /P  in  Eq . (C-5)  then  represents  the  power 

ol  io 

back-off  at  0-dB  normalized  input. 

4.  Initialization 

The  program  is  initiated  by  specifying  as  inputs  the  satellite  termi- 
nal characteristics  listed  in  Table  C-l  and  assigning  specific  frequencies 
and  powers  (ERP)  to  the  n links,  either  directly  or  from  an  initialization 
plan.  The  bandwidth  W of  the  repeater  is  defined  as  the  frequency  differ- 
ence between  the  highest  and  the  lowest  assigned  frequency  links  plus  half 
the  sum  of  the  data  rates  of  these  two  links.  The  initial  frequency  plan 
for  the  remaining  carriers  is  then  a proportionate  spacing  based  on  the 
data  rate  of  each  link.  Fixing  the  end  points  of  the  repeater  bandwidth 
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reduces  the  number  of  variables  in  the  frequency  optimization  by  two.  The 
frequency  variable  used  in  the  optimization  procedure  is  the  difference 


r 


between  successive  frequencies,  i.e.,  Af  f - f , Af  = f - f . and 

1 1 o 2 2 1 


The  ERP  are  assigned  in  proportion  to  the  ratio  of  the  data  rate  and 
recei vi ng-termina 1-ga i n- to-t empera t lire  ratio  (G/T)  such  that  the  sum  of 


the  power  levels  equals  a constant.  Thus,  if  there  are  n links,  the  ERP 


th 


for  the  i link  is  determined  as 


R /(G/T) 
i i 


s 


S 


t i 11  to 

R /(G/T)k 

i.  ' 


(C-8 ) 


where  S is  the  constant  representing  the  equivalent  power  of  a single 
to 


carrier.  The  program  is  set  up  so  that  any  desired  value  can  be  assigneci 


to  S.  . 
to 


A subroutine,  IPERM,  is  included  in  SYSCON  to  run  the  program  with 
the  different  possible  orderings  of  the  links  in  the  repeater. 


Example 


Suppose  we  want  to  determine  the  value  of  to  yield  a given  value, 


say  1/p,  for  P.  /P  in  Figure  A-7.  Thus,  we  require 

J in  omax  " ’ ' 


p • G • G 

st  twt 


who  re 

P 

omax 

is  the  maximum  output 

power  of 

the  limitei 

r i n 

Figure  C-l.  Th 

power 

at  the 

input  to  the  limiter 

P a nd 

in 

the  t ransmi tier 

ERP  S are  re- 

t o 

lat  ed 

by  the 

range  equation 

G 

• G • 

2 

X • 1, 

P = S 
in  to 


sr  a 


( 1 rid  ) 
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Thus,  S 


is  given  by: 


to 


s 

to 


P 

max 


G • G • G * G 
sr  a twt  st 


or,  in  decibels, 


S (dBW) 
to 


P -G  -G  -G  -G-L 
max  sr  st  twt  a 


p + 20  log 


For  an  up-link  frequency  of  8.0  GHz  and  a synchronous  satellite  with  the 
following  repeater  characteristics 


G 

sr 

16.8 

dB 

Satellite  receiving  antenna  gain 

G 

st 

17.0 

dB 

Satellite  transmitting  antenna  gain 

P 

max 

13.0 

dBW 

Satellite  EHP 

G 

a 

69.0 

dB 

Gain  of  the  preamplifier  preceding 

the  limiter 

L 

1.0 

dB 

Miscellaneous  link  losses 

G 

twt 

27.0 

dB 

Gain  of  the  TWT  amplifier 

we  obtain 

S =74  dBW 

to 

if  p = 10  dB, 

i . e . , P . 

/P 

equal  to  -10  dB  is  required. 

in  omax 


5.  Expression  for  the  Error  Probability 
a.  General 

The  general  expression  for  the  bit-error  probability  at  the 
output  "f  the  i*  channel  is  given  bv  [Appendix  B,  Kq.  ( B-20 )] 


( C-9 ) 


7 K 


who  re 


2 

r = 
i 


y\2  s .w  s 

/ j i k i k / j i j i j 


(C-10) 


k = 1 
k 1 


j = l 


t h 

represents  the  equivalent  SNR  or  2E/n  at  the  output  of  the  i channel. 

o 

2 2 *2 

The  SNRs  i , . , and  p at  the  receiver  can  be  determined  by 

i i k 1 j 

considering  the  down-link  transmission  of  the  signals: 


2 

2P  • (o  /T  \ • A2 
r \ ri  ri  / i 

2P  • (g  /T  \ • A2 
r \ ri  ri/  i 

i 

p ■ Ig  n ,il 

K R 

r \ ri  ril 

i 

K 

1 + — T 

R 

16  s P 

i 

t J 

/ \ 2 

/ \ 

2P  • (G  n ) ■ A, 

2P  • (G  ,/T  • A 

2 

r \ ri  ri  / k 

r \ ri  ri/ 

ik 


1 + — T 
16  s 


P • (C  /T 
r t ri  ri 


K R 


(C-l 1) 


(C-12) 


2P  • (G  /T  ) • A2  2P  • (g  /r  ) • A2 

r \ ri  ri/  m n p r \ ri  ri/  m n 


ij 


1 + — T 
16  s 


P • G 'T 
r l ri  ri 


K R 


(C-13 ) 


and 


p = 
t 


i = 1 


n 1/2 


M-.' 


j + 


i = l 


i = 1 


(C-M) 
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I 


where 


P = satellite  ERP  referenced  to  the  ground  terminal 
r 

th 

G /T  = figure  of  merit  of  the  i receiver 
ri  ri 

K = Boltzman's  constant 

T - satellite  noise  termperature  in  degrees  Kelvin 
s 

th 

s = power  of  the  i signal  at  the  repeater  input 

1 

2 

Y~  = limiter  parameter 
i 

R = bit  rate  of  the  quadriphase  signal 
i 

th 

A = amplitude  of  the  i signal  at  the  limiter  output 
i 

th 

A = amplitude  of  the  k signal  at  the  limiter  output 
k 

A = amplitude  of  the  third-order  cross  product  gener- 

m n p 

ated  by  the  mixing  of  signals  on  links  m,  n,  and  p. 
Here  the  A coefficients  are  identical  with  the  h coefficients  discussed 


in  Appendix  A.  They  are  defined  as: 


a 


( C - 1 5 ) 


where  a.  is  the  amplitude  of  the  i signal  at  the  limiter  input,  and  - 
is  the  noise  power  at  the  input  to  the  limiter.  The  approximations  in 
Eqs.  ( C — 1 1 ) through  (C-13)  refer  to  the  fact  that  in  FDMA  the  retransmitted 
up-link  noise  is  usually  negligible  compared  to  the  down-link  noise. 


The  satellite  KRP  P referenced  to  the  ground  terminal  and  the  i 
r 

signal  power  s at  the  input  to  the  satellite  repeater  are  determined  as 
i 


fo 1 lows : 


p = p 
r max 


u (1 


(■I'd) 


(C-16) 


h. 


2 

s = s • G d d 
i ti  sr  — 

(1  -d  ) 


P = satellite  ERP 
max 


(C-  17) 


} = wavelength  of  the  midband  frequency  for  down-link 


t ran sm i ss ion 


1,  = miscellaneous  losses  (down-link) 

d 

d = distance  between  the  ground  station  and  the  satellite 

G = satellite  receiving  antenna  gain 
sr 

th  . 

s = ERP  of  the  l ground  terminal, 
ti 

By  combining  F.qs.  (B-38)  through  (B-48)  in  Appendix  11,  the  fol- 
lowing expressions  are  obtained  for  S and  S^.  in  Eq . (C-10): 


(*-JqL)('-J^)Hv)' 


T - min  T , for  ? - i,k 
min 


(C-1K) 


S = — 
ij  T, 


T = min  T.,  for  i1  = i,m,n,p  , (C-19) 
min  " 


r 


where 
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<3 
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J 
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T 

i 

2 

R. 

= bit  interva 

1 of 
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th 

i 

1 i nk 

R 

i 

= bit 
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th 

i 

1 i nk 

“i 

= radian  frequency 

of 

the  i 

th 

1 i nk 

= radian  frequency 

of 

tho  k 

tb 

link 

= radian  frequency 

of 

t he  .j 

th 

cross 

.1 


Tn  Eqs.  (C-18)  and  (C-19)  the  indices  i,  m,  n,  and  p indicate 
the  links  involved,  while  the  index  ,j  is  an  index  on  the  i ntermodulation 
product  generated  by  signals  from  links  m,  n,  and  p. 


Three  subroutines  are  used  in  SYSCON  to  calculate  the  error  rate 
by  Eq . (C-9):  AMP,  SPECOLP,  and  ANORM.  Subroutine  AMP  calculates  the 
amplitudes  of  the  signals  and  the  cross  products  at  the  output  of  the  band- 


pass limiter;  SPECOLP  determines  the  power  contribution  at  the  output  of 
th 


the  i receiver  from  all  the  other  FbMA  signals  (adjacent-channel  inter- 
ference) and  third-order  cross  products;  and,  finally,  ANORM  calculates 
the  bit-error  rate  for  each  channel  by  using  Eq.  (C-9). 


b.  Amplitudes  of  Limiter  Output  Signals  and  Intermodulation 
Components 


Program  AMP  calculates  the  amplitudes  of  tho  signals  and  inter- 
modulation product;-  by  Eq . (C-15).  The  infinite  integral  is  approximated 
by  a finite  integral  in  the  interval  from  0 to  U,  where  U is  the  upper 
limit  of  integration  and  is  determined  such  that  the  truncation  error  in 
the  integral  is  negligible.  Equation  (C-15)  is  first  normalized  through 
substitution  of  t lie  following 


8 2 


1 


2 2 

1 +Yi 


2 2 
- + V 


c = 
2 


c = 
3 


where  c = min  C > . This  yields 


1 

A , . . . = - J 

P,  P„  P 2 
12  n ; i = 1 i 

o 


p (a.  • z/  c)  exp (-  -f  /) 

+ exp(“T?2  t) 


(C-20 ) 


Since 


Hj;  ' 


£ 1 


i = 1 i 


and  at  z = 1 


-z  -16 

e e -8 

= = 2.85  X 10 

4 4 


the  integral  converges  to  a limit  very  fast  as  z increases.  It  was  found 
that,  under  all  circumstances,  the  contribution  of  the  integral  in  the 
region  z = 4 to  oo  is  typically  of  the  order  ol  0.0001  percent.  So  the 
infinite  integral  in  Eq.  (C-20)  can  be  approximated  by  integration  in  the 
region  z = 0 to  4.  Thus,  =4  in  Eq.  (C-20). 

The  Bessel  functions  J and  J in  Eq . (C-20)  can  be  calculated 

o 1 


from  their  polynomial  approximations 


1 3 
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x <1  3 


-:s  ^ 

2 , 4 

J ( x ) = 1 - 2.24999  97  (x/3)  + 1.26562  08(x/3) 

o 

6 8 
- 0.31638  66  (x/3)  + 0.04444  79  (x/3) 


10 

- 0.00394  44 (x/3) 


, 12 

0.00021  00 (x/3) 


| e|  ■'  5 X 10 


(C-2  1) 


3 < x < “ 

-1/2  -1/2 
J (x)  = X r cos  0 Y (x)  = X i sin  9 

o o o o O o 

2 

f = 0.79788  456  - 0.00000  077(3/x)  - 0.00552  740(3/x)“ 
o 

- 0.00009  512(3/x)'{  + 0.00137  237(3/x)1 

5 . 6 

- 0.00072  805 (3/x)  + 0.00014  476 (3/x)  + ,~ 


e < 1.6  y 10 


(r-22) 
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x - 0.78539  816  - 0.04166  397(3/x) 


0.00003  95  1 (3/x)2  + 0.00262  573(3. 'x)3 


0.00054  125(3/ x)  - 0.00029  333(3/x)° 


0.00013  558 (3/x) 


6 


I e|  " 7 y 10 


-8 


-3  x 3 


- 1 1 2 
x Jx(x)  = - - 0.56249  985 ( x/3 ) + 0.21093  573(x/3) 


- 0.03954  289 (x/3 ) ^ + 0.00443  319(x/3)8 


10  1 9 

- 0.00031  761 (x/3)  + 0.00001  109 (x/3)  “ 


1.3  x 10 


-8 


, „ -1/2 
J ( X ) = X f cos 
1 1 1 


3 i x 


fj  = 0.79788  456  + 0.00000  156(3/x)  + 0.01659  667(3/x)' 


0.00017  105(3/x)'i  - 0.00249  511(3/x)  ' 


0.00113  653 ( 3/x ) ° - 0.00020  0.33 (3/x) 6 


? <4  y 10 


-8 


(C  -23 ) 


(C-24) 


(0-25) 
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6 

0.0029  166  (3/x)  + c 


e < 9 y 10 


-8 


(C-26) 


and  is  calculated  by  the  recurrence  relationship 


J (z)  + J (z)  = — J (z) 

V - 1 +1  Z V 


(C-27 ) 


The  finite  integral  of  Eq . (C-20)  is  first  divided  into  either  two  or  four 

subintervals,  depending  upon  the  maximum  value  of  a //c;  then  each  sub- 

i 

interval  is  evaluated  by  the  16-point  Gaussian  quadrature  method. 

The  value  of  a / rc  in  J and  J influences  the  accuracy  with 
i o 1 

which  the  integral  in  Eq.  (C-20)  can  be  calculated  by  the  Gaussian  quad- 
rature method.  A larger  value  of  a /J  c means  J and  J will  have  more 

i o 1 

zero  crossings  in  the  interval  (0,4).  Hence  a polynomial  of  higher  degree 
is  needed  to  approximate 


i = 1 

and  hence  more  points  will  be  needed  in  the  Gaussian  quadrature  method. 

Subroutine  AMP  takes  this  into  consideration.  For  a Tc  1,  the  integral 

i 

interval  (0,4)  is  divided  into  two  subintervals,  (0,2)  and  (2,4).  For 
a.y/c  > 1,  the  interval  (0,4)  is  divided  into  four  subintervals,  (0,0.5), 
(0.5,1),  (1,2),  and  (2,4).  Each  subinterval  is  then  evaluated  by  the 
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Subroutine  AMP  evaluates  Eq . (C-20)  with  the  aid  of  three  sub- 

rou tines: 


FUNCTION  AJO(X) — To  calculate  Bessel  function  J . 

<» 

FUNCTION  AJI(X) — To  calculate  Bessel  function  .J^. 

SUBROUTINE  QG 16 — To  perform  16-point  Gaussian  quad- 
rature integration. 

Subroutine  AMP  has  been  thoroughly  tested.  Values  of  the  inte- 
gral calculated  from  this  program  were  compared  with  results  from  an 
earlier  SRI  program  that  utilizes  the  Newton-Cotes  integration  formula. 

The  results  from  these  two  programs  were  found  to  be  in  excellent  agree- 
ment (to  at  least  five  figures).  The  program  using  the  Gaussian  quad- 
rature method  runs  very  much  faster  than  the  one  using  the  Newton-Cotes 
formu la . 

c.  Frequencies  of  the  Cross-Product  Components 

Subroutine  CALXSM  is  used  in  SYSCON  to  calculate  the  frequencies 
of  all  third-order  intermodulation  products.  Note  that  there  are  two  types 
of  third-order  cross  products.  These  are  generated  at  frequencies 


is  produced  as  a result  oi 

f , respectively. 

.1 
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2 f.  - f.  = f , 
l .j  21 


and 


f - f 

.1  k 


1 1 1 


The  cross  product  at  the  frequency  2f  - f 

i .j 

mixing  of  two  signals  at  frequencies  f,  and 


The  cross 


d.  Calculation  of  and  Sjj 

Subroutine  SPECOLP  calculates  the  integrals  in  Eqs.  (C-18)  and 
(C-19).  Note  that  Eq.  (C-18)  is  a subset  of  Eq . (C-19),  involving  only 
two  products  instead  of  four.  By  multiplying  out  the  four  products  in 
Eq.  (C-18),  S can  be  expressed  as  a sum  of  four  integrals 


T 

/•min 


T • S 
1 U 


/ 1 I \ J*  m 1 n 

r T - (r + r + f + f ,cos K >T 

\ m n p i / J 


1 1 1 

+ 


T T T T T T 
( m n m p mi 


T T T T T T 
n p n i pi 


2 

T COR 


s/aiu  . t|c 


-( 


T T T T T T T T T T T T 
mnp  mni  npi  mpi 


3 

t cos  / Au) . ' \d 


T T T T 
m n p i 


t co  s 


K )■ 


T - min  T,,  for  l - i,m,n,p  . (C-28) 

m j n i K 


Each  term  can  be  integrated  analytically, 
used  in  programming  Eq . (C-28), 
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Table  C-2  gives  the  expressions 


INTEGRALS  USED  IN  THE  SUBROUTINE  SPECOLP 


here  is  the  same  as  £j.  in  Eq.  (C-28). 


dices.  As  Aui . approaches  zero,  the  integrals  simplify,  with  the  cosine 


J 

terms  approaching  one.  For  values  of  Auu  T close  to  zero,  the  expres- 

.1  min 

sions  shown  in  Column  3 of  Table  C-2  are  used. 

To  find  the  expression  for  S , it  is  necessary  only  to  set  1/T 

ik  p 

and  1/T  equal  to  zero  in  Eq . (C-28). 

n 

Equation  (C-28)  is  used  as  follows  for  the  various  types  of 
third-order  cross  products: 

• If  the  cross  product  under  consideration  is  generated 

as  a result  of  mixing  of  signals  other  than  the  desired 

signal,  Eq . (C-28)  is  valid  for  determining  the  spectrum 

overlap  resulting  from  the  cross  product  of  the  type 

f + f - f . For  the  cross  product  of  the  type  2f  - f , 
m n p m n 

we  need  only  to  set  1/T  equal  to  zero  in  Eq . (C-28). 

P 

• If  the  cross  product  under  consideration  is  generated  as 
a result  of  mixing  of  signals  involving  the  desired  sig- 
nal, then,  for  the  type  f + f - f cross  product,  S 

m n p i ,j 

is  determined  by  setting  either  1/T  , 1/T  , or  1/T  equal 

pm  n 

to  zero  (depending  upon  which  of  the  three  indices  m,  n,  p 

is  equal  to  i).  For  the  cross  product  of  the  type 

2f  - f , we  set  1/T  equal  to  zero,  as  explained  above, 
m n p 

and  also  either  1/T  or  1/T  equal  to  zero  (depending 
m n 

upon  which  of  the  two  indices  m or  n is  equal  to  i). 

This  procedure  is  necessary,  as  explained  in  connection 
with  Eq.  (B— 47),  when  considering  the  effects  of  cross 
products  at  the  receiver  that  are  produced  because  of 
mixing  of  signals  with  the  desired  signal. 


90 


Approach  for  Selecting  Power  and  Frequency 


6 . 

The  approach  taken  for  optimization  in  SYSCON  is  to  optimize  the 

norm  N or  N first  with  respect  to  either  the  frequency  or  the  power, 
a b 

by  holding  one  of  them  fixed,  and  then  to  repeat  the  procedure  on  the 

parameter  that  was  held  fixed.  The  reason  for  doing  this  is  that  fre- 

2 2 

quency  and  power  are  easily  separable  in  Eq . (C-10).  The  SN'Ils  p , p , 

i i k 

-2  - 
and  o are  functions  of  the  transmitter  ERPs,  while  S and  S are 
i j ik  l.j 

dependent  on  their  frequencies. 


The  program  is  initiated  by  first  assigning  specific  frequencies 
and  powers  to  the  n links,  either  directly  or  from  the  initialization 
plan  discussed  in  Section  4.  The  probability  of  error  for  each  of  the 
n links  is  then  calculated.  The  next  step  is  to  minimize  either  of  the 
two  norms,  which  can  be  selected  at  the  start  of  the  program. 

The  problem  of  minimizing  the  norm  is  equivalent  to  finding  the 
minimum  of  a function  of  n variables.  A standard  technique  for  solving 
a problem  of  this  type  is  the  method  of  steepest  descent.  It  consists  of 
picking  a point,  calculating  the  partial  derivatives  at  this  point,  and 
proceeding  in  the  direction  of  the  steepest  slope  by  a specified  amount. 
After  the  step  has  been  taken,  the  function  is  evaluated  to  see  whether 
it  is  less  than  its  previous  value.  If  it  is  less,  the  procedure  is  con- 
tinued; if  not,  a smaller  step  is  taken.  If  a minimum  exists,  it  will  be 
found  e ve  n t ua 1 1 y • 

The  algorithm  used  in  SYSCON  for  both  power  and  frequency  optimiza- 
tion is  an  improved  version  of  the  steepest  descent  method. 

a.  Description  of  the  Optimization  Approach 

Methods  of  minimization  of  a function  of  several  variables  can 
be  divided  into  two  groups:  direct  search  methods,  which  do  not  depend 
explicitly  on  evaluation  or  estimation  of  the  gradient  vector  of  the 


St  l 


objective  function,  and  gradient  methods,  which  require  derivatives. 
Since  no  derivatives  are  required  by  the  direct  search  methods,  they  are 
applicable  to  functions  with  discontinuous  derivatives.  If  derivatives 
can  be  evaluated,  the  gradient  methods  are  usually  far  more  efficient. 

In  the  satellite  FDMA  problem,  analytic  expressions  for  the  derivatives 
of  the  norm  exist;  therefore,  only  gradient  methods  for  minimization  of 
the  norm  were  considered. 

In  the  discussions  below,  U(c)  is  a function  of  k variables, 

[UC^)]  is  the  minimum  value  of  U(c),  and  i stands  for  the  it!' 
min 

di rect ion . 


1)  The  Newton  Method 


If  we  expand  U(c  + Ac)  in  a Taylor  series  about  c,  we  have 


U(?  + AS)  = u(  C)  + 711  (3)a3  + ... 


(C-29) 


where  U is  the  gradient  vector.  Differentiating  the  Taylor  series  again, 
we  obtain 


•7U(C  + A^)  = 7U  ( c ) + iiA2  + ... 


(c-:u)) 


where  H,  the  Hessian  matrix,  is  given  by 


K) 


(C-31 ) 


For  c + A?  to  be  the  minimizing  point,  C , '7U(3>  + Ac)  would  be  0,  so, 

mi  n 

neglecting  higher-order  terms, 


— U 

AC  = - H 0) 


(C-:i2) 


9 2 


If  we  are  dealing  with  a quadratic  function,  this  incremental  change  will 
take  us  to  the  minimum  in  only  one  iteration. 


r 


When  ('  is  not  quadratic,  one  can  try  the  iterative  scheme 


j + 1 


cj  - rf  V 


(C-.Ti) 


where  H is  the  inverse  of  the  Hessian  matrix  at  the  j iteration.  This 

scheme  has,  however,  several  disadvantages.  H must  be  positive  definite; 

1~  j 

otherwise  divergence  could  occur.  In  particular,  -H  U might  not  point 
downhill.  Also,  H may  be  locally  singular;  if  so,  its  inverse  does  not 
exist.  This  scheme  was  tried  for  the  optimization  of  the  norm.  The 
computations  of  H,  the  second  derivative  matrix,  and  its  inverse  were  very 
time  consuming. 


2)  Steepest  Descent 

Instead  of  Eq . (C-33),  one  can  try  the  iterative  scheme 


y.l 


1 + 1 


a s 


(C-.'i  1) 


where  s'  is  (hopefully)  a "downhill''  direction  of  search,  and  a'1  0 is 

the  step  size  or  a scale  factor  chosen  to  minimize  U (f'1  + cf1*1). 

The  most  obvious  choice  for  s1  is  the  steepest-descent 
direction  at  J1.  From  Fq . (C- 29),  a first-order  change  in  the  objee  • i 
function  is  given  by 

~ T „ 

AU  = -ru  AJ 

I f A-1  = ,’s1,  where  o'  > 0 is  fixed  and 

the  s'  minimizing  “l  or  maximizing  -AF  i-  . > • n 


D-A047  044 


UNCLASSIFIED 


STANFORD  RESEARCH  INST  MENLO  PARK  CALIF  F/G  17/2.1 
INVESTIGATION  OF  THE  USE  OF  FREQUENCY-DIVISION  MULTIPLE  ACCESS  — ETC(U) 
DEC  73  P C JAIN*  V E HATFIELD*  J K LEUNG  DCA100-72-C-0033 

NL 


2<*2 

ADA047044 


END 

DATE 

FILMED 

I -78 

DOC 


where  ‘ 1 is  the  Euclidean  norm  defined  as 


1/2 


The  s'5  in  Eq.  (C-36)  is  the  negative  of  the  normalized  gradient  vector. 
Although  - 7U/!|^U!l  provides  the  greatest  local  change,  success  of  the 
steepest  descent  method  is  highly  dependent  on  the  scaling  factor  cy"1 . 

3)  Algorithm  Used  in  SYSCON  for  FDMA  Optimization 

The  selected  algorithm  is  a modified  version  of  the  steep- 
est descent  method.  The  objective  function  U corresponds  to  the  norm  N, 
and  ® is  either  the  frequency  or  the  input  signal  amplitudes.  Equation 
(C-36)  of  the  steepest  descent  method  is  modified  as  follows: 

~ j ~ j - 1 

_ i tu  vu 

SJ  = (1  - 3)  + 3 — , (C-37 ) 

l|7UJl|  l|^J_1|| 

to  locate  the  direction  of  the  minimum.  Comparison  of  this  s'1  with 
Eq . (C-36)  shows  that  the  gradient  selected  is  the  weighted  sum  of  the 
two  gradients.  For  3 = 0 we  obtain  the  same  result  as  with  the  steepest 
descent  method. 

The  step  size  a [in  Eq.  (C-S-l)]  is  first  set  to  be  propor- 
tional to  | ® ^ | , and  is  further  adjusted  automatically  until  a minimum  is 
located  along  the  s^  direction;  then 


sJ  + 1 


= e 


aJSJ 


and 

u(«J  + 1)  < u(cj) 


(C-38) 


is  satisfied. 

The  iteration  is  repeated  until  some  prescribed  stopping 
criterion  is  satisfied.  The  stopping  criteria  consist  of: 

• The  maximum  number  of  allowable  iterations  has 
been  reached. 

• Either  one  of  two  inequalities  has  been  satisfied, 

- u(»J)l  * <J 

- *hl  = e_  , 

& 

where  e and  e have  prescribed  values. 

1 2 

A simplified  flow  chart  of  the  optimization  Subroutine  STEEP  is  shown  in 
Figure  C-3  and  is  discussed  in  Volume  One,  Section  IV-E. 


7 . Derivatives  of  the  Norm 

For  power  and  frequency  optimization,  we  need  the  derivatives  of  the 

norm  with  respect  to  the  signal  amplitudes  and  frequency  differences.* 

Thus,  if  norm  N is  used,  we  need  to  determine 
a 


As  mentioned  in  Section  4,  the  variable  used  in  frequency  optimization 
is  the  difference  between  successive  frequencies,  i.e.,  Afj  = f ^ - f , 

Af 2 = f2  ~ fj,  and  so  on,  instead  of  absolute  frequencies  fJ(  f 2 , and  so 
on.  This  procedure  reduces  the  number  of  variables  in  the  optimization 
procedure  by  two. 
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Here  a and  Af  represent,  respectively,  the  desired  amplitude  and  frequency 


differences. 


a.  Determination  of  9Pej/9a 

— 

It  can  easily  be  shown  from  Eq.  (C-9)  that 


9 P 


ei  e 


-p./2(l+h)  p 


9a 


/2t7 


pi 


(1  t h) 


9h 


3/2  3a 


3p, 


<11  h) 


1/2  3a 


(C-41 ) 


where 


h = P2wS  ■ u + P2  ■ ■ 
Z—J  ik  lk  ij  ij 

k=  1 

k^i 


(C-42 ) 


and 


n 3p 

9h  „ ik 

9a  ik  9a 

k=l 
k^i 


9p. 


+ 2 > S 


ij 


i j 9a 


(C-43) 


The  derivatives  of  the  SNRs  with  respect  to  signal  amplitude 

can  be  expressed  in  terms  of  the  corresponding  derivatives  of  the  limiter 

output  amplitudes,  by  Eqs.  (C-ll)  through  (C-13) . Thus,  in  general,  we 

would  require  the  derivatives  3A  /9a , 3A  / 9a , and  9A  /'a,  which, 

i k m n p 
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in  turn,  can  be  calculated  from  Eq.  (C-15).  Since  derivatives  of  Bessel 
functions  can  be  expressed  in  terms  of  other  Bessel  functions,  they  lead 
to  integrals  similar  to  Eq . (C-15)  and  can  be  evaluated  by  using  the  ap- 
proach discussed  in  the  calculation  of  Eq.  (C-20) . In  SYSCON  the  deriva- 
tive of  the  norm  N or  N with  respect  to  the  signal  amplitudes  is  calcu- 
a b 

lated  by  evaluating  the  analytic  expressions  in  Subroutine  DNORMA. 


b.  Determination  of  AP„.j/AAf 


£1_ 


From  Eq . (C-9)  we  can  derive  the  following  expression: 


AP 


ei  1 _ e 
2 


-p  /2(  1+h) 
i 


AAf 


/2n 


(1+h) 


Ah 


3/2  AAf 


where  h is  defined  by  Eq.  (C-43)  and 


. 11  AS  AS 

Ah  _ 2 ik  ^ ,2  ij 

AAf  / -<  1 ik  AAf  AAf 


k=l 

kri 


The  derivatives  of  Eq . (C-45)  can  also  be  expressed  as 


AS 


AS 


ik 


and 


AAf  A Au). 


AS  As 

ij  = iJ 

AAf  AAuu  . 


A Aid 


AAf 


A Au) 


AAf 


(C-44 ) 


(C-45) 


(C-46) 


(C-47 ) 


Analytic  expression  for  the  derivatives  in  Eqs.  (C-46)  and  (C-47) 

can  be  developed  from  Eqs.  (C-18)  and  (C-19).  In  SYSCON,  the  derivatives 

of  the  norm  N or  N with  respect  to  the  frequency  differences  are  calcu- 
a b 

lated  in  Subroutine  DNORMF . 
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8.  Summary 

A simplified  conceptual  flow  chart  diagram  of  SYSCON  is  shown  in 
Figure  C-4 . The  required  input  data  are  identified.  The  satellite  and 
terminal  parameters  are  those  given  in  Table  C-l.  At  the  start  of  the 
program  it  is  essential  to  provide  the  initial  ordering  of  the  links  in 
the  satellite  repeater,  as  well  as  selection  of  the  norm  and  the  indica- 
tion whether  optimization  of  the  norm  is  desired  Also  necessary  is  the 
information  whether  other  possible  combinations  of  the  links  are  to  be 
considered. 

If  the  limiter  parameters  Yj  and  y are  not  given  in  the  input  data, 
they  are  determined  in  Subroutine  PINOUT.  If  they  are  known,  this  pro- 
cedure is  bypassed. 

Next,  the  program  determines  the  initial  starting  values  for  the 
ERPs  and  the  up-link  frequencies  of  the  transmitters.  Three  options  are 
available.  First,  the  values  and  frequencies  can  be  provided  in  the 
input  data.  Second,  they  can  be  determined  from  the  initialization  plan 
discussed  in  Section  4.  A third  option  allows  the  possibility  of  pro- 
viding only  the  values  for  the  transmitter  ERPs  at  the  input;  the  initial 
frequency  assignment  is  then  assumed  to  be  uniform  within  the  satellite 
repeater. 

The  program  then  calculates  the  probability  of  error  P for  each 

ei 

link  by  using  Subroutines  AMP,  SPECOI.P,  and  ANORM  and  also  determines  the 
value  of  the  selected  norm.  If  minimization  of  the  norm  was  indicated  at 
the  input,  it  is  performed  next  in  Subroutine  STEEP.  The  order  of  opti- 
mization (i.e.,  power  first  and  then  frequency,  or  vice  versa)  and  the 
number  of  optimizations  must  be  provided  at  the  input.  Subroutine  STEEP 
(Figure  C-3)  is  used  for  both  power  and  frequency  optimization.  The 
values  for  the  stopping  criteria  are  already  included  in  the  program,  but 
they  can  easily  be  changed  if  desired.  The  values  used  in  the  program  are 


i 
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INPUT  DATA 


SATELLITE  PARAMETERS 
TERMINAL  PARAMETERS 
ORDERING  OF  THE  LINKS 
SELECTION  OF  THE  NORM 


PERMUTATION  ON  THE  ORDERING 
OPTIMIZATION  OF  THE  NORM 


S A - 1975-40 


FIGURE  C-4  SIMPLIFIED  FLOW  CHART  OF  SYSCON 
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Number  of  iterations  M = 10 


ex  = 10 


e2  = 0.001  x 


i = l 


At  the  end  of  each  optimization  (power  or  frequency)  the  values  of  P . , 

ei 

S . , f , and  the  norm  are  printed, 
ti  l 

After  the  last  optimization,  the  program  can  be  exercised  again  with 
the  next  combination  of  the  links,  if  that  is  desired.  The  possible  com- 
binations are  determined  in  Subroutine  IPERM. 

The  numerical  results  obtained  by  using  SYSCON  are  discussed  in 
Volume  One. 


v. 

\ 
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(1-dB  back-off)  than  is  customary  in  practice. 

To  obtain  the  expression  for  the  error  rate  at  the  input  of  the  FDMA  links,  it 
was  necessary  to  derive  general  analytic  expressions  for  the  limiter  output  signal , 
the  intermodulation,  and  the  noise  components,  when  n signals  are  transmitted 
simultaneously  through  the  satellite  repeater.  The  expression  for  the  bit  error  rate 
was  then  derived  by  assuming  digital  quadriphase  modulation  of  the  FDMA  carriers  and 
taking  into  consideration  the  presence  of  other  FDMA  carriers  causing  adjacent-channel 
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13.  Abstract  (Concluded) 


A general  comparison  of  the  three  major  multiple-access  alternatives — 
FDMA,  TDMA,  and  SSMA--with  respect  to  selected  performance  criteria  that 
are  particularly  important  for  the  military  environment  and  for  operation 
with  a mix  of  users  indicated  that  FDMA  performs  much  better  than  past 
analyses  had  shown.  With  the  same  satellite  power  and  RF  bandwidth,  FDMA 
was  found  to  offer  nearly  as  much  satellite  throughput  as  TDMA  and  con- 
siderably more  than  SSMA. 

The  report  consists  of  three  volumes.  Volume  One  provides  a summary 
of  the  study  which  includes  a description  of  the  analysis  approach,  docu- 
mentation of  the  pertinent  equations,  numerical  results,  and  conclusions. 
Detailed  analysis  and  investigation  of  the  problem  areas,  as  well  as 
derivations  of  analytical  expressions,  are  contained  in  the  appendices  in 
Volume  Two.  Volume  Three  provides  a technical  assessment  of  the  suitability 
of  PSK/FDMA  for  operation  with  a mix  of  users  and  compares  its  performance 
with  other  multiple-access  alternatives,  in  particular,  with  tl  t of 
time-division  multiple  access  (TDMA)  and  spread-spectrum  multiple  access 
(SSMA)  . 


